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Resumo 


Obtivemos famílias de soluções estáticas e esfericamente simétricas para buracos 
negros numa brana com vácuo, constante cosmológica não nula e imersa num bulk 
5-dimensional assintoticamente anti-de Sitter. Para o caso das geometrias assintotica- 
mente de Sitter, o único membro da família obtida foi a geometria Schwarzschild-de 
Sitter extrema. Para o caso de geometrias assintoticamente anti-de Sitter, obtivemos 
toda uma família de soluções, onde cada elemento desta família é determinado por uma 
constante C. Estudamos o comportamento destas soluções próximo e longe do hori- 
zonte de eventos. A seguir, analisamos a evolução de um campo escalar não massivo 
no exterior dos buracos negros obtidos. Verificamos analiticamente e numericamente 
o comportamento do potencial efetivo em função da coordenada tartaruga próximo ao 
horizonte de eventos. Para qualquer geometria ou solução desta família observamos que 
o campo escalar não massivo decai com o tempo, e o tipo de decaimento depende do 


valor da constante C. Assim, toda a familia de soluções obtida mostrou-se estável. 


ill 


Abstract 


We found classes of static and spherically symmetric black hole solutions on the 
brane for a vacuum configuration, non-null cosmological constant, embedded into a 5- 
dimensional asymptotically anti-de Sitter bulk. In the asymptotically de Sitter case, 
the only solution obtained was the extreme Schwarzschild-de Sitter geometry. In the 
asymptotically anti-de Sitter case, a class of solutions was obtained where a particular 
geometry is fixed by a constant, C. We studied the behaviour of these solutions near and 
far from the horizon. Next, we studied the time evolution of a scalar field outside the 
black hole. For the effective potential as function of tortoise coordinate, the analytical 
results near the horizon agree with numerical calculations. The scalar field decays with 
time and the particular form of decay depends of the value for C. Thus, all solutions 


in the obtained class are stable. 
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Capitulo 1 


Introducao 


As palavras teoria e relatividade, quando juntas, referem-se a uma das idéias mais 
populares da ciéncia, embora pouco compreendida pela maioria das pessoas. Quando 
Albert Einstein publicou sua Teoria da Relatividade em artigos nos anos de 1905 e 1916, 
nao poderia prever o impacto de suas idéias na visão de mundo e na vida das pessoas, 
muito menos a popularidade de sua teoria. A Teoria da Relatividade modificou as 
noções de espaço e tempo, tornando-os interdependentes, de tal forma que não podemos 
tratá-los isoladamente, mas como um espaço-tempo. Ambos os conceitos, a partir de 
então, não são absolutos como na Mecânica Newtoniana, mas relativos, dependem do 


observador. 


Outra grande contribuição de Einstein foi relacionar geometria com matéria ou 
energia, isto é, matéria determina a geometria e vice-versa. No entanto, sua equação, 
proposta em 1916, foi modificada quando a constante cosmológica, Asp, foi adicionada 
para que sua teoria estivesse de acordo com a idéia de um universo estático em grande 
escala. Entretanto, a constante cosmológica foi abandonada quando Edwin Hubble 
constatou a expansão do universo na década de vinte do século passado e, somente 
no final do século XX, foi ressuscitada para tentar esclarecer a expansão acelerada do 


universo. 


Em 1916, logo depois da publicação da Teoria da Relatividade Geral, um físico 
alemão chamado Karl Schwarzschild obteve a solução exata das equações de Einstein 
no vácuo com constante cosmológica nula para uma geometria esfericamente simétrica. 
Tal solução ficou conhecida como solução de Schwarzschild. Os elementos da métrica 


de Schwarzschild, no sistema de coordenadas esférico, divergem para dois valores da 
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coordenada radial, r: em r = 0 e em r = 2m. A última divergência pode ser removida, 
desde que uma mudança para um sistema de coordenadas apropriado seja feita. Como 
veremos no capítulo 2, um sistema de coordenadas poderá não ser eficaz para descrever 
completamente todo um espaço-tempo. A singularidade em r = 0 é conhecida como 
singularidade física ou essencial e o valor 2m corresponde ao raio do horizonte de even- 
tos, que é uma superficie de mao única, ou seja, depois de penetrar nesta superficie não 
há como sair dela, nem mesmo a luz consegue. Por isso, esta solução descreve o buraco 


negro de Schwarzschild. 


Como foi dito acima, a Teoria da Relatividade tranformou dois entes distintos num 
único só. O espaço-tempo quadridimensional é onde e quando eventos ocorrem. Entre- 
tanto, no mesmo século, nos anos vinte, Kaluza e Klein nos trabalhos [1, 2] propuseram 
a existência de uma dimensão extra compacta na tentativa, sugerida pelo primeiro, de 
unificar o eletromagnetismo com a gravitação. Podemos pensar em dimensões com- 
pactas através do conhecido exemplo da mangueira no jardim. Vista do alto, uma 
mangueira parece uma linha, um objeto unidimensional, entretanto se chegarmos perto, 
veremos que este objeto além de comprimento, possui largura e espessura, ou seja, é 
um objeto tridimensional. Nos anos 90, a teoria que buscava, e ainda busca, unificar as 
forças da Física é a Teoria de Cordas, através do uso de 6 dimensões espaciais extras 
compactas (veja por exemplo [3]). Dessa forma, temos um espaço-tempo com 9 dimen- 
sões espaciais e 1 dimensão temporal, este espaço-tempo 10-dimensional é chamado 
de bulk. Atualmente existe a possibilidade, desde que modelos teóricos mostrem-se 
realistas, de detecção indireta de dimensões extras através de colisões de partículas 


elementares no LHC (Large Hadronic Collider), como pode ser visto em [4]. 


A unificação da Fisica torna-se necessária pois em altos níveis de energia efeitos 
quânticos não são desprezíveis, o que conduz à unificação das duas teorias mais bem 
sucedidas do século XX, a Teoria da Relatividade e a Mecânica Quântica. Além disso, a 
Teoria de Cordas é uma candidata a resolver o problema da hierarquia, que é a diferença 
entre a escala de massa eletrofaca, My ~ 10°GeV, e a escala de massa de Planck 4- 
dimensional, M, ~ 10!°GeV, ou seja, é a diferença entre as forças eletrofraca e a força 
gravitacional, sendo M, = (era = onde Gy é a constante gravitacional Newtoniana. 
Obter uma teoria de campos que envolva o modelo padrão e a gravitação seria difícil 
devido a essa diferença de escalas. Uma opção para resolver o problema da hierarquia 
é o uso de dimensões extras. Neste contexto, modelos de mundo brana surgem como 
tentativa de resolver este problema com auxílio de dimensões extras. O modelo tipo 


Randall-Sundrum II, que será usado neste trabalho, consiste de uma brana quadrimen- 
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sional imersa num bulk anti-de Sitter 5-dimensional com simetria Z2. Todos os campos, 


exceto o gravitacional, estão confinados na brana. 


Bronnikov et al. mostra em [5| como obter uma ampla classe de soluções que 
descrevem buracos negros estáticos esfericamente simétricos na brana. A partir da 
função A(r) = gu da métrica pode-se obter uma família de funções B(r) = —1/grr 
e, logo em seguida, fixando um valor para uma determinada constante de integração 
C a solução fica completa. Cada membro da família corresponde a um valor de C 
e a uma determinada geometria. Neste trabalho, utilizamos os A(r)’s das geometrias 
Schwarzschild-de Sitter e Schwarzschild anti-de Sitter para obter uma família de funções 
B(r) numa brana com vácuo e constante cosmológica não nula. Embora para o primeiro 
caso a única solução obtida tenha sido a já conhecida solução Schwarzschild-de Sitter 
extrema, para o segundo caso, obtivemos numericamente uma família de soluções na 
brana. Estudamos as características desta família longe e próxima do horizonte de 


eventos, determinando, por exemplo, a gravidade superficial. 


O passo seguinte foi verificar a estabilidade da família de soluções obtida. Esta- 
bilidade das soluções das equações de Einstein é objeto de estudo desde a década de 
50 do século passado quando Regge e Wheeler, no trabalho [30], estudaram a solução 
de Schwarzschild submetida a uma perturbação. Ambos verificaram que a solução de 
Schwarzschild é sempre estável, ou seja, a perturbação diminui com o tempo. A família 
de soluções obtida, a partir de um A(r) Schwarzschild anti-de Sitter, mostrou-se estável 
para o caso de uma perturbação escalar para quaisquer valores de C. Todo o cálculo 
de estabilidade, que envolveu a evolução do campo escalar ®, foi realizado com auxílio 


de cálculo numérico, por meio do problema das coordenadas duplas nulas. 


Podemos, então, descrever todos os passos contidos neste trabalho da seguinte forma: 
o capítulo 2 refere-se à geometria diferencial, ou seja, todo instrumental matemático 
necessário para se ter alguma compreensão da Teoria de Relatividade Geral; o capítulo 
3 refere-se aos mundos branas, com destaques aos modelos Randall-Sundrum I e II e 
à projeção das equações de Einstein na brana; o capítulo 4 refere-se a buracos negros 
estáticos com simetria esférica, com destaques às definições de horizontes de eventos 
e Killing e a construção de soluções desse tipo; o capitulo 5 refere-se às familias de 
soluções de buracos negros estáticos e esfericamente simétricos na brana obtidas a partir 
do A(r) Schwarzschild-de Sitter e Schwarzschild anti-de Sitter; o capítulo 6 refere-se à 
estabilidade da família de soluções obtida no capítulo 5 quando submetida a um campo 


escalar não massivo; o capítulo 7 refere-se às conclusões finais deste trabalho. 


Capitulo 2 
Geometria e Relatividade Geral 


Teoria Geométrica da Gravitação poderia ser o nome da Teoria da Relatividade Geral 
tamanha a importância da Geometria Diferencial na sua formulação. Por isso, neste 
capítulo um pouco de geometria será exposto com intuito de nos fornecer uma base 
para compreensão dos capítulos subseqiientes. Conceitos como variedade, métrica entre 


outros serão apresentados. 


2.1 Introduzindo a geometria 


Um conceito fundamental em geometria diferencial é o de variedade. Tal conceito 


captura a idéia de um espaço que pode ser curvado e ter uma topologia complicada, 


porém localmente se assemelha ao R”. Para se ter um uma definição mais rigorosa de 
variedade é necessário ter em mente outros conceitos tais como mapa, difeomorfismo, 


sistema de coordenadas, atlas e espaços tangentes. 


Um mapa, à, entre dois conjuntos M e N é a relação que assegura para cada elemento 
de M um elemento em N. Tal relação pode ser expressa por 9: M — N. Dado dois 
mapas ¢: A— B e y : B > C, podemos definir a composição Pod: A > C pela 
operação (bo q)(a) = Y(¢(a)). Dessa forma, a € A, (a) E Be (Wod)(a) E C. Existem 
mapas do tipo injetor, sobrejetor e bijetor. Um mapa ¢ é injetor se cada elemento de 
N tem no máximo um elemento de M correspondente. Se o mapa for sobrejetor, cada 
elemento de N tem no mínimo um elemento de M correspondente. Caso o mapa seja 
injetor e sobrejetor, ao mesmo tempo, ele se chamará bijetor. Neste caso, podemos 


definir o mapa inverso 9 !:N — M por (¢-'o da) = a. O conjunto M é chamado 


CAPITULO 2. GEOMETRIA E RELATIVIDADE GERAL 5 


de dominio do mapa q e o conjunto de pontos em N, que é mapeado por M, é chamado 


imagem de q. 


Uma bola aberta, na topologia induzida pela métrica usual, é o conjunto de todos 


pontos 7 no R” tal que |z — y|< r, para algum y € R” fixo e r € R, onde |z| é uma 


norma usual. Um conjunto aberto no R” é construído de uma arbitrária união de bolas 


abertas, em outras palavras, V C R” é aberto se, para qualquer y € V, há uma bola 


aberta centrada em y que está completamente dentro de V. Se um mapa entre conjuntos 


abertos do R” é p vezes diferenciável ele é dito C? (está sendo usada a definição usual de 


diferenciabilidade entre abertos do R”). Caso um mapa seja infinitamente diferenciável, 


será dito liso e indicado por CS. Dado duas variedades, M e N, e os mapas dg: M > N 
egt: N — M lisos e invertíveis, M e N serão chamados difeomórficos e ¢ será chamado 


de difeomorfismo. 


Um sistema de coordenadas consiste de um sub-conjunto U de um conjunto M, que 


está relacionado com o R” através de um mapa injetor ¢: U — R”, tal que a imagem 


ọ(U) é um aberto no R”. Um atlas é uma coleção de sistemas de coordenadas {(Ua, 


da) | que satisfaz duas condições: 


1. A união de Us é igual a M, isto é, Ua cobre todo M. 


2. Se dois sistemas de coordenadas se interseccionam, Ua N Ug 4 Ø, logo o mapa 
(da © 5°) leva os pontos de ¢g (Us N Ug C R”) em da (Ua N Ug CR?) e todos 


esses mapas devem ser lisos onde são definidos. 


A importância em se falar de atlas e sistema de coordenadas reside no fato de 
que muitas variedades não podem ser inteiramente cobertas por um único sistema de 
coordenadas. Como exemplo, podemos citar o espaço-tempo de Schwarzschild, que não 
é bem descrito pelas coordenadas (t,r,0, 4) em uma região que inclui o horizonte de 


eventos. 


Introduziremos, logo em seguida, estruturas características de uma variedade: os 
espaços tangente e cotangente. O espaço tangente V, num ponto P de uma variedade 
M, é formado por operadores que levam funções na variedade a números reais. Isto é, 


dado um operador derivada direcional d/dA, onde À é um parâmetro ao longo de uma 
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curva em M, sua aplicação a uma função f definida na variedade pode ser escrita como 


df — dat 


ENE ax ul, (2.1) 


onde ð, = 0/02” é a derivada parcial. A base desse espaço é formada por operadores 
derivada parcial, 0,,, que são bases de um espaço vetorial de derivadas direcionais que 
é identificado como o espaço tangente (veja [15] capítulo 2). A partir de (2.1), vemos 


que vetores contravariantes dx"/d\ = X” são elementos desse espaço. 


O espaço cotangente ou dual V*, num ponto P da variedade M, é formado por 


operadores que levam vetores do espaço tangente a números reais, ou seja, w : V > R. 


Tais operadores são 1-formas (covetores) que podem ser indicadas pelo gradiente df de 
uma função f 

d df 

df(—) = =, 2.2 
o) = (2.2) 
onde d/dA é um vetor tangente a uma curva de parâmetro À na variedade e pode ser 


escrito com auxílio de (2.1) como 


d da! 

DE qd (2.3) 
e, dessa forma, (2.2) fica 

df =O, fdr": (2.4) 


Podemos, então, interpretar o gradiente das coordenadas, dx", como a base do espaço 
cotangente, pois, por definição, uma base de um espaço tangente, ê,, e uma base de um 
espaço cotangente, à”, devem satisfazer ú!ê, = ô}. Esta relação é satisfeita quando ó, 


e dx! são bases dos espaços tangente e cotangente, respectivamente. Ou seja, 
CL Oy = ô$. (2.5) 


Uma 1-forma pode, por conseqiiéncia, ser escrita como w = wyda" e, além disso, pode- 


mos observar que vetores covariantes w, são elementos desse espaço dual. 


Em uma variedade diferenciável M, tensores são definidos de forma natural como 


aplicações multilineares de vetores e covetores em um ponto de M até R. O conjunto 


de aplicações de ordem (k,l) forma o espaço tangente de tensores 


T: V* 8- GV OVO--OV SR (2.6) 
k l 
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com k índices covariantes e | índices contravariantes. Considerando o conjunto de 
espaços vetoriais de um dado tipo em cada ponto da variedade, obtém-se um campo 
vetorial. É importante ressaltar que sem estruturas adicionais não existe relação entre 
espaços tangentes em pontos distintos. Entre estas estruturas adicionais está derivada 


covariante, que será discutida na seção 2.4. 


Qualquer tensor covariante simétrico de ordem 2 define uma métrica g : V@V 5 R, 
ou em termos de suas componentes em um sistema de coordenadas, g(x). Uma 
variedade equipada com uma métrica é chamada variedade Riemanniana. A métrica 
fornece uma noção de causalidade e quadricomprimento e é usual representá-la através 


da distância infinitesimal, ds, entre dois pontos x” e x!” + da”, 
ds? = g (x)dx”dz”. (2.7) 


A equação (2.7) é conhecida como elemento de linha. O quadrado do comprimento ou 


norma de um vetor X” é definido como 
X? = gula) X!X”, (2.8) 


onde vetores do tipo espaço têm X? < 0, vetores do tipo tempo têm X? > 0 e vetores 


do tipo luz ou nulo têm X? = 0. 


Podemos representar a métrica na sua forma canônica e, dessa forma, suas compo- 


nentes são indicadas por 
Guy = diag(+1, +1, ..., —1, —1, 0, ..., 0), (2.9) 


onde diag significa os elementos da diagonal principal da matriz que representa gv. Se 
representarmos o número de +1’s por t e o número de -1’s por s, a diferença entre s — t 
nos fornece a assinatura da métrica. Quando uma métrica tem s = 0 ela é chamada 
de Euclidiana ou Riemanniana e quando tem t = 1 ela é dita Lorentzina. Quando diag 
possui zeros, a métrica será degenerada. Na Relatividade Geral, métricas de interesse 


são do tipo Lorentziana e não degeneradas. 


Uma solução das equações de Einstein, uma métrica, pode não ser bem comportada 
numa determinada região do espaço-tempo quando escrita num sistema de coordenadas 
qualquer. Devido ao fato de que uma variedade pode não ser coberta inteiramente por 
um único sistema de coordenadas, uma mudança para um outro sistema de coordenadas 


pode ser realizada e, dessa forma, a métrica poderá descrever a região anteriormente não 
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descrita de forma adequada. ‘Tal procedimento chama-se extensão analítica da métrica. 


2.2 Simetrias em Relatividade Geral 


Entre os conceitos mais importantes em geometria diferencial estão os conceitos de 
pullback e pushforward. Tais conceitos nos levarão a definições de isometria e campo de 


Killing e, a partir daí, às noções de espaços estáticos e estacionários na seção seguinte. 


Vamos considerar duas variedades M e N, que podem não ter o mesmo número de 


dimensões, com sistemas de coordenadas x” e y“, respectivamente. O mapa ¢: M > N 


e a função f : N > R podem ser representados de forma ilustrativa através da figura 
2.1: 


Mo N a 
e ina! 
E E 


R” R 


Figura 2.1: Pullback de f por ¢ (figura baseada em [15] capítulo 5) 


Podemos compor à com f e construir o mapa (fo ¢): M — R, uma função de M em 


R. Tal operação pode ser indicada por 


pf =(f © 4), (2.10) 


e tem o nome de pullback de f por q. Como o nome sugere, “puxamos” a função f 
para trás, de N para M. Caso tivéssemos uma função em M, g: M — R, então 


não poderíamos fazer o pullback de g por ¢, não poderíamos compor à com g. Isto é 
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percebido através do sentido das setas. Dessa forma, não podemos “empurrar” a função 


g para frente, de M para N. 


Como vetores podem ser pensados como operadores em espaços vetorias tangentes 
em cada ponto da variedade e mapeiam funções lisas (funções contínuas e infinitamente 
diferenciáveis) em números reais, podemos definir o pushforward de um vetor. Sendo 
V(p) um vetor num ponto p em M, o pushforward de V por ¢ em um ponto ¢(p) em 


N será 


(OVI) = VoF). (2.11) 


Então o vetor “empurrado” para frente atuando numa função f em N será igual a esse 


mesmo vetor atuando numa função “puxada” para trás em M. 


Bases para um espaço vetorial tangente num dado ponto das variedades M e N 
são dadas por operadores derivadas parciais, como já foi mencionado. Num ponto 
da variedade M, temos um espaço vetorial cuja base será 0, + e num ponto da 


variedade N temos outro espaço vetorial cuja base será ða = 5. Podemos relacionar 


— Oxo 


as componentes de V = V”, em M com as componentes de (9*V)ºda em N aplicando 


a operação de pushforward a uma função lisa com auxílio da regra da cadeia 


(DV) df V" On sf), 
= V"Ə (f od), 
Oy* 
Ve Oe 2.12 
0a f (2.12) 
Escrita desta forma, a operação de pushforward lembra uma operação de mudança de 
coordenadas, pois em (2.12) dy" /dx” faz o papel de uma matriz de mudança de coorde- 
nadas. Entretanto, tal operação será uma generalização para mudança de coordenadas 


apenas quando M e N são a mesma variedade. 


Não podemos usar a operação de pullback para vetores, podemos usá-la para 1- 
forma, que são vetores duais. Como uma 1-forma é um mapa linear de vetores a 
números reais, o pullback de uma 1-forma w em N, w, pode ser definido por sua ação 
em um vetor V em M igual à ação de w num vetor onde a operação de pushforward foi 


realizada 
(hw) (V) = w(9*V). 


Novamente, usando a regra da cadeia, a operação de pullback em formas lembra uma 


operação de mudança de coordenadas. 
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Um tensor Tá,...o, covariante em N é uma aplicação linear do produto direto de | 


covetores ou 1-formas em R. Assim, podemos usar a operação de pullback não somente 


em formas mas também em tensores covariantes arbitrários. Então, a definição da 
operação de pullback em tensores covariantes é sua ação em vetores que sofreram a 


operação de pushforward 
(¢.T)(VY,V®,..,V) = TE VO VO, VD). (2.13) 


Da mesma forma, um tensor contravariante S“!“* pode sofrer a operação de pushfor- 


ward 


(ES (WD w, ...,0) = Slow, bw, ..., dw). (2.14) 


As representações das matrizes das operações pullback e pushforward para tensores de 


ordem elevada podem ser escritas como 


Oy Oy™ 
CDE) u.p = pm om OA (2.15) 
oy oð 
* QI. Qk __ y é y“ 1- Hk 
(¢*S) Spa GRE (2.16) 


para um tensor T de ordem (0,/) e um tensor S de ordem (k,0). Tensores mistos 
podem ser “puxados” e “empurrados”, desde que à seja um difeomorfismo, isto é, seja 
liso e possua uma inversa, $!, igualmente lisa. Como um exemplo, para um tensor 


Trote, mem M, a operação de pushforward para pode ser representada por 


(PT)(w™, ut VOLVO = T(d.0™, pw P, THV ®, (VOS), 
(2.17) 
onde ws são 1-formas em N e os Vs são vetores em N. Podemos representá-la em 
componentes da seguinte forma 


QUO coro OD a cna, 
Prt Data Oech dy dy na (2.18) 


e aa 
O pullback de um tensor misto Tits, |, em N pode ser representado por 


(TD, cpu, VO, VO) = TP Ww, 2 OW, HV, HV), 
(2.19) 


onde os ws são 1-formas em M e os V:s são vetores em M. Podemos representá-la 
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em componentes da seguinte forma 


dr” = Ox” dy Oy 
[o ao __ HiH 
Eia eaa a a (2.20) 


Como já foi dito, se representarmos o pullback e o pushforward em componentes, 
obteremos matrizes de transformação de coordenadas quando M e N são a mesma va- 
riedade. Quando q é um difeomorfismo e as variedades M e N são a mesma, obtemos 
um tipo de transformação de coordenadas que é chamada de transformação “ativa”. 
A transformação de coordenadas tradicional é chamada de “passiva”. Vamos conside- 


rar, como exemplo, uma variedade n — dimensional M com funções de coordenadas 


at: M > R”. Podemos mudar de coordenadas introduzindo novas funções de co- 


ordenadas y“ : M — R”, mantendo a variedade fixa. Outra forma de realizar essa 


mudança de coordenadas é usar um difeomorfismo à : M > M obtendo as novas co- 


ordenadas através de pullbacks e pushforwards (¢.x2)" : M > R”, movendo os pontos 


na variedade e depois calculando suas coordenadas. 


Um difeomorfismo permite-nos fazer operações de pullback e pushforward em ten- 
sores. Podemos, dessa forma, comparar valores de um campo tensorial em diferentes 
pontos na variedade como, por exemplo, um campo Tt», |, (x) calculado num ponto 
p e num ponto ¢(p), onde a operação de pullback foi aplicada, [TAr #* 1,..1,(@(p))]- 
Assim, um outro tipo de derivada é obtido, uma derivada que mostra como campos 
tensorias mudam sobre um determinado difeomorfismo. Entretanto, precisamos não 


apenas de um difeomorfismo, mas de uma família de difeomorfismos de um parâmetro, 


ġı. Tal família pode ser pensada como um mapa liso R x M > M, de tal forma que 


para cada t € R, à; é um difeomorfismo e Ps o dy = Ps4+, sendo que ġo é a identidade. 


Famílias de difeomorfismo de um parâmetro podem ser pensadas como originárias 
de campos vetoriais e vice-versa. Em cada ponto na variedade, sob uma família de 
Qı, temos curvas. Essas curvas cobrem toda variedade e, dessa forma, podemos definir 
campos vetorias, V(x), que são o conjunto de vetores tangentes a essas curvas em 
cada ponto, calculados em t = 0. De forma contrária, podemos definir uma famila de 
difeomorfimos de um parâmetro através de campos vetoriais. Dado um campo vetorial, 
V(x), podemos definir as curvas integrais, x(t), como sendo as curvas que resolvem 
a equação 

—— = Va). (2.21) 


Soluções para (2.21) em variedades são garantidas por teoremas que podem ser vistos, 


por exemplo, em [18] capítulo 4. 
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Podemos obter uma famíla de difeomorfismos parametrizada por t, a partir de um 
campo vetorial V”(x), e calcular a rapidez com que tensores variam conforme percorrem 
curvas integrais, que são soluções de (2.21). Para um t qualquer, podemos definir a 


variação de um tensor por 


PAN Ride Hk V1... (p) = Pts [oe et OUD) Se Vi.. (p) ; (2.22) 


Podemos representar esta variação através da figura 2.2: 


ESET ae 
ee aaa TI «(P)] 


PTC ED + 
ne 


Figura 2.2: Variação de um tensor Tite, , (figura baseada em [15] capítulo 5) 


Dessa forma, podemos definir a derivada de um tensor ao longo de um campo vetorial. 


Tal operação chama-se derivada de Lie 


A, Pra ibe es 
Ly THM yg = Hi (SE at) (2.23) 
A derivada de Lie é independente do sistema de coordenadas, é uma mapa de um tensor 
do tipo (k,l) a um tensor do tipo (k,l), é linear (V (T + S) = VT + VS, onde Te S 
são tensores) e obedece a regra de Leibniz (V(T 8 S) = (VT) 8 S +T @(VS)). 


Para facilitar a discussão sobre a derivada de Lie e outras operações relativas à ela, 

vamos adotar um sistema de coordenadas x“ em que x! é um parâmetro ao longo das 
E E ; . eee.) . 

curvas integrais. Assim, o campo vetorial toma a forma V = 55, ou seja, 0 campo 

tem as componentes V” = (1,0,0,...,0) e um difeomorfismo de uma família de um 

parâmetro, ¢;, produz a transformação de coordenadas de x“ a y” = (al +t, x?,..., a”). 


Conseqtientemente, a matriz de pullback fica 


ae (2.24) 
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e as componentes do tensor que sofreu a operação de pullback de ¢:(p) a p ficam 
Dia Pa (Co) dg a a (2.25) 


Nesse sistema de coordenadas, a derivada de Lie de um campo tensorial é expressa por 


o 
Ly TEk i E agit a (2.26) 
ou para um campo vetorial U” (x) 
OU" 


A derivada de Lie pode também ser expressa por um comutador 
Ly JU" = [V,U], (2.28) 
sendo este comutador dado por 
MUP SOUS — U” 8, V”. (2.29) 


Vamos fazer o uso de (2.28), com o auxílio de (2.29), para o caso da métrica guv- 


Obtemos, então, 


V°Voguw + (VuV gx + (VV gu 
= VV tV Va 
= aval 2.30) 


Ly Gi 


onde V, é a derivada covariante de guv- 
Dizemos que um difeomorfismo, ¢, é uma simetria de algum tensor, T, se depois 
que aplicamos a operação de pullback neste tensor ele permanece o mesmo, ou seja, 


oT = T. (2.31) 


Se uma famíla de simetrias de um parâmetro, œ+, é gerada por um campo vetorial V“(z), 
por (2.27) usando (2.31) temos 
Err EU: (2.32) 


Uma consequência desta simetria é o fato de que o tensor T é simétrico sob uma família 
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de difeomorfismos de um parâmetro e podemos sempre obter um sistema de coordenadas 


onde as componentes de T são sempre independentes de uma coordenada 


As componentes de T não dependem de xt. 


Uma importante simetria é aquela relacionada à métrica, ou seja, O.Juv = Juv- 
Um difeomorfismo que gera tal simetria é chamado de isometria. Se uma família de 
isometrias de um parâmetro é gerada por um campo vetorial V”(x), este campo é do 


tipo Killing. A condição que V” seja um campo de Killing é a seguinte: 
Lv uv = 0, (2.34) 


ou usando (2.30) 
ViVi =0. (2.35) 


A equação (2.35) é conhecida como equação de Killing. Então, se um espaço-tempo 
tem um campo de Killing, podemos escrever a métrica que o representa independente 


de uma coordenada. 


2.3 Espaços estacionários, estáticos e esfericamente 


simétricos 


Para espaços com dimensão n > 2, podemos ter mais campos de Killing do que 
dimensoes. Isto se deve ao fato que um conjunto de campos de Killing pode ser line- 
armente independente, mesmo que num dado ponto da variedade os vetores de Killing 
sejam linearmente dependentes, e, dessa forma, pode gerar um outro campo de Killing 


(desde que a combinação linear seja a coeficientes constantes). 


Um espaço de dimensão n tem simetria máxima quando tem o número de campos de 


Killing linearmente independentes associados a ele dado por: n(n + 1)/2. Por exemplo, 


no R” temos duas classes de isometrias de um parâmetro: translações e rotações. Há 


n translações, uma para cada direção, e n(n — 1) rotações (para cada direção há n — 1 
direções que podemos rodar, e para evitar uma contagem repetida dividimos por 2, 


ou seja, fazer uma rotação em torno de z no sentido horário e anti-horário é a mesma 
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coisa). Assim temos 
—1 1 
n4 mn a nn ) (2.36) 


campos de Killing no R”. Como exemplo de geometrias com simetria máxima temos o 


espaço-tempo de Minkowski, cuja métrica é dada por 
ds?’ = dt? — dx” = dy = dz’, (2.37) 


onde t é a coordenada temporal e (x,y,z) são coordenadas cartesianas. Outros casos 


importantes são o espaço-tempo de Sitter (dS), cuja métrica é dada por 


2 2\ 71 


a a? 


e o espaço-tempo anti-de Sitter (AdS), cuja métrica é 
r2 r2 =] 
ds? — (1 + Z) dt? — (1 + =) dr? — r2(do? + sen7dd’), (2.39) 
a a 


sendo t a coordenada temporal, r a coordenada radial, 0 e ¢ são as coordenadas angu- 


lares. A constante a é positiva e está relacionada com a constante cosmológica por 


3 
a? = +, 2.40 
TE (2.40) 


onde Asp > O no caso de Sitter e Ayp < 0 no caso anti-de Sitter. 


Um espaço-tempo que possui simetria esférica têm campos de Killing linearmente 


independentes que satisfazem as relações de comutação abaixo: 


VD, V7] = Ver 


[V®, V3] =- V® 


ve ve) =e), (2.41) 


Como exemplo, podemos citar o espaço-tempo de Schwarzschild, que possui a métrica 
2 2m 2 2m\ 2 2/1092 20742 
ds* = | 1 — — (dt — | 1 — — dr“ — r*(d0" + sen*0d¢*), (2.42) 
r r 


onde m é uma constante real positiva, que pode ser interpretada como a massa. 
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Um espaço-tempo é dito estácionário se ele admite um campo de Killing do tipo 
tempo [14]. Portanto, uma métrica é estacionária se existe um sistema de coordenadas 


onde ela é independente do tempo, isto é, 


O Ju 
=), 2.43 
a70 (2.43) 
sendo x° é uma coordenada do tipo tempo. Se definirmos um campo vetorial X” = ô}, 


logo a derivada de Lie da métrica será 
Lx Qu = X°VoGu + (ViX*) gn + (VA) guar; 


LxIw = 05 Vouw i 0, (2.44) 


por (2.43). Como Lyg,, é um tensor, então ele será zero em qualquer sistema de 
coordenadas e X” é um campo de Killing, de acordo com (2.34), do tipo tempo. De 
forma inversa, dado um campo de Killing do tipo tempo X”, sempre existe um sistema 


de coordenadas que satisfaz (2.44). 


Um espaço-tempo é dito estático se ele pode ser folheado por uma família de hiper- 
superfícies ortogonais a um campo de Killing do tipo tempo, de acordo com [14]. Num 
espaço-tempo estático, há um sistema de coordenadas adaptado ao campo de Killing 
do tipo tempo de tal forma que a métrica é independente da coordenada temporal e 
nao há termos cruzados, isto é, goa = O. Portanto, para falarmos de espaço-tempo es- 
tático, precisamos do conceito de campo vetorial ortogonal a uma hipersuperfície (uma 
hipersuperfície em uma variedade de dimensão n é simplesmente uma subvariedade de 


dimensão n — 1). A equação de uma família de hipersuperfícies é dada por 
E) =a, (2.45) 


onde diferentes valores de a correspondem a diferentes hipersuperfícies. Dado dois pon- 
tos P e Q com coordenadas x” e (x! + dx"), respectivamente, na mesma hipersuperfície 
S, por (2.45) podemos escrever em primeira ordem de aproximação em torno do ponto 


P 
of 


Fa! + dx") = f(a") + TT 


da”. (2.46) 


Entretanto, f(x” + dx!) = f(x") pois os pontos P e Q estão na mesma hipersuperfície. 


Então, o segundo termo de (2.46) é nulo no ponto P. 
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Se definirmos um campo vetorial covariante 7, como 


Of 


Tu Deo 


(2.47) 


logo 
Hat Qua): de = 0 (2.48) 


no ponto P e nos diz que n, é ortogonal à dx”. Consequentemente, n, é ortogonal à 
hipersuperfície S pois dx” está em S e é chamado de campo vetorial normal à S em P. 


Qualquer outro campo vetorial X” ortogonal à S será proporcional à 7,,, isto é, 
AP SAH, (2.49) 


sendo A(x) um fator que depende do ponto onde calculamos. Pode-se mostrar que 
A(x) = X? = XºX, (veja [14] capítulo 14) e assim (2.49) fica 


of 
— 2 
X =X 


(2.50) 

Uma métrica estática não contém termos cruzados que envolvem a coordenada do 
tipo tempo, isto é, uma solução estática não contém termos do tipo dxzľdx”. Para 
entender isso, considere o intervalo entre os eventos (x°, x!, °, x3) e (x° + dx?, xt + 


num istem rden ; rm 
da er dado sistema de coordenadas. Dessa forma, ds? fica 


ds? = goo(dx°)? + 2goidx°dz' + g1 (dx')?. (2.51) 
Sob uma reversão temporal 1º > x’? = —2°, (2.51) fica 
ds? = goo(dx°)? — 2g da" da! + gn(dr!)?. (2.52) 


Como numa solução estática ds? é invariante sob uma reversão temporal, go, deve ser 
zero para que (2.51) e (2.52) coincidam. Da mesma forma, se tivéssemos usado x? + dx? 


ou z? + dxê, go» ou gogz deveriam se anular igualmente. 


Dado um espaço-tempo estacionário e um campo de Killing do tipo tempo, isto é, 
X" = 6h, logo, 
Xp = Gi = Gio = Jop (2.53) 


X? = X, X” = Jou% = goo- (2.54) 
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Então (2.50) fornece com auxílio de (2.53) e (2.54) 


of 
Jon = 00 5 > (2.55) 


para algum campo escalar f. Quando u = 0, (2.55) fica Of /Ox° = 1 e depois de 
integrá-la obtemos 
f=2°+h(2), 


onde h é somente uma função das coordenadas do tipo espaço e o indice a vai de 1 a 3. 


Agora, vamos fazer uma transformação de coordenadas 
z? —> z” =2°+4 A(z’), L rl =r. (2.56) 


Como no caso estacionário, a métrica não depende da coordenada do tipo tempo 


1 
Juy 
Ox’? 


=0. (2.57) 


Nesse sistema de coordenadas, X'” = 6g e com o auxílio de (2.55) e (2.56) obtemos 


On. 
Jou = oo gr (2.58) 
Para o caso = 0 temos 
Goo = goo. (2.59) 
e para o caso u + 0 conseguimos 
Joa = Uy (2.60) 


isto é, não há termos cruzados do tipo dx°dx* e a solução, dessa forma, é estática. 


2.4 Conexão, curvatura e equações de Einstein 


Nas seções anteriores foram introduzidas uma variedade diferenciável e uma métrica 
Lorentziana como estruturas básicas para descrever o espaço-tempo. Na variedade são 
definidas de maneira natural funções, curvas parametrizadas, espaços tangentes, vetores 


e tensores. 


Pontos distintos em M não compartilham, em geral, do mesmo espaço vetorial 


tangente, exceto o R”, que tem para cada ponto um espaço tangente associado a ele 
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igual ao próprio R”. Dessa forma, o operador derivação do R” não pode ser o operador 
de derivação de uma variedade qualquer. Para uma variedade qualquer, queremos 
um operador V, que seja independe do sistemas de coordenadas, seja linear em seus 
parâmetros, obedeça a regra do produto de Leibniz e se transforme como um tensor. 
O operador derivada covariante aplicado a um tensor contravariante, Xº, numa dada 
direção y, é dado por 

Van OP TRA, (2.61) 


e para o caso de um tensor covariante, Xa, é dado por 


Via = 0X = Di Xp (2.62) 
Se exigirmos que V, Xº seja um tensor do tipo (1,1), T3, deve se transformar de acordo 


com i 
a OE Ox? LT y OL? OF? On 


bro dao Ax'® ðr °° Ax'B Ax" Ox Ax?’ 


Se o segundo termo do lado direito de (2.63) fosse zero, teríamos a usual lei de trans- 


(2.63) 


formação de tensores do tipo (1,2). Portanto, a presença do segundo termo do lado 
direito revela que a lei de transformação é inomogénea, ou seja, l'3, nao é um tensor. 
Qualquer objeto que se transforme de acordo com (2.63) é chamado de conexão afim. 
Uma variedade com uma conexão afim contínua é dita variedade afim. A equação (2.61) 


num outro sistema de coordenadas pode ser escrita como 


Ve SOX Dn dg (2.64) 
sendo 
DX = O (52º) 
ô ta 
E (ae) 
T oe eee ass 
: ox’? 
X’? = re (2.66) 


O segundo termo de (2.65) faz com que Ox nao seja um tensor. Substituindo (2.65), 
(2.63) e (2.66) na equação (2.64), com uma troca apropriada de índices mudos veremos 


que o segundo termo de (2.63) é cancelado com o segundo termo (2.65), fazendo com 
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que o operador derivada covariante transforme-se como um tensor. 


Uma conexão em que a derivada covariante da métrica é zero é chamada de “com- 
patível com a métrica”. Existem infinitas conexões deste tipo. Porém, existe uma única 
conexão compatível com a métrica e com torção nula. Esta é a conexão usada na 


relatividade geral, denominada conexão de Christoffel 
a 1 ad 
By = 59° (p95 + A958 — A596) - (2.67) 


Introduzida a conexao, podemos definir a curvatura. Para a conexao de Christoffel, 
o tensor curvatura ou tensor de Riemann é dado por 
R” as = 0,85 — O51 + Da oy — Dá e (2.68) 


Com o uso de contrações do tensor de curvatura, define-se várias quantidades impor- 


tantes: o tensor de Ricci 


Rw ae g” Rie: (2.69) 


e o escalar de curvatura ou escalar de Ricci 


R= GG Ria (2.70) 


A introdução de matéria é fundamental na descrição da gravidade. Na relatividade 
geral, isso é feito (em geral) através de campos definidos na variedade, associados a um 
tensor de energia-momento, Tau. Estes campos obedecem a equações de movimento 
que em geral podem ser expressas em termos de relações tensoriais. Nestas equações, 
as derivadas espaciais são derivadas covariantes com relação à conexão definida pela 


métrica. 


Em 1916, Einstein propôs como equação de movimento para as componentes da 
métrica a relação 

Gm = ie Su = = (2.71) 

A equação tensorial acima, um sistema de equações não lineares para as 10 componentes 

independentes da métrica, é conhecida como equação de Einstein e foi modificada pelo 

próprio Einstein quando a constante cosmológica, Asp, foi adicionada para que sua 


teoria estivesse de acordo com a idéia de um universo estático em grande escala. Então, 
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a equação tensorial tornou-se 


8TG 
Gw Eu AsDIuv = a Tur (2.72) 


A constante 87Gy /c* é escolhida de forma que no limite não relativístico sejam obtidos 
os resultados previstos pelas leis de Newton. É bastante comum na literatura, como 
por exemplo em [25], o uso das chamadas “unidades geométricas”, em que Gy = c=1. 


Neste sistema de unidades, usado deste ponto em diante, o fator 8rG y/c torna-se 87. 


Capitulo 3 
Universos Brana 


Neste capitulo, discutiremos possíveis dimensões extras, exploraremos modelos de mun- 
dos brana, com destaque aos modelos Randall-Sundrum, e, na seção final, obteremos a 


projeção das equações de Einstein generalizadas na brana. 


3.1 Dimensões extras 


A utilização de dimensões extras para descrever fenômenos físicos foi feita pioneira- 
mente por Kaluza e Klein nos trabalhos [1, 2] e, atualmente, é feita pela Teoria de 
Cordas. Busca-se, conseqiientemente, a influência destas dimensões extras no nosso 


mundo 4-dimensional e modificações em Teorias Físicas. 


Segundo a lei da Gravitação Universal de Newton, o módulo da força entre duas 
partículas de massas my e Mz, separadas pela distância r, é dado por Gunmima/r?, 
onde Gy é a constante gravitacional newtoniana. O fato da força depender somente 
da distância relativa entre os corpos pode ser creditado como uma conseqiiéncia da 
isotropia do espaço, ou seja, as leis físicas são independentes da direção, são as mesmas 
em qualquer direção. Quando adota-se h/27 = c = 1 (onde h é a constante de Planck ec 
é a velocidade da luz), Gy = (10"*cm)? é inversamente proporcional a MZ, onde M, = 
10'°Gev é a massa de Planck. Dessa forma, temos uma grandeza (Gy = 10em?) 


muito pequena, incapaz de ser observada atualmente. 


Para o caso de objetos astronômicos, objetos macroscópicos, a lei da gravitação 
diminui com 1/r? e funciona muito bem. Quando são feitas medidas em pequenas dis- 


tâncias, a influência de outras forças, como a força elétrica por exemplo, interfere na 


22 
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obtenção de resultados. Isso é devido ao fato de que a lei da gravitação é fraca a curtas 
distâncias. Os melhores dados obtidos por [6] mostram que abaixo de 0, Imm não pode- 
mos discernir entre a força gravitacional e as demais forças. Então, a lei 1/r? funciona 
bem entre escalas astronômicas, distâncias da ordem do sistema solar, e 0, Imm. Caso 
dimensões extras existam, sua importância seria notada em medidas abaixo de 0, Imm. 
Além disso, com uma dimensão espacial extra, a gravitação seria notada em todas di- 
mensões e diminuiria com 1/rê. De uma forma geral, num mundo n-dimensional, a 
lei da gravitação diminuiria com 1/r”7? e seria inversamente proporcional à massa de 
Planck n-dimensional. Assim, a massa de Planck 4-dimensional seria relacionada com 
a n-dimensional por 


M$ = MV, (3.1) 
onde V é o volume das n — 4 dimensões adicionais. 


Em 1921, Kaluza [1| na sua tentativa de unificar o eletromagnetismo com a gravita- 
ção propôs um espaço-tempo com 4 dimensões espaciais, isto é, 3 dimensões infinitas e 
uma quarta compacta presa num círculo de tamanho re na dimensão extra. Klein num 
trabalho em 1926 [2], determinou que o tamanho desse círculo seria da ordem da escala 
de Planck, 10% cm, algo totalmente não observável. Dessa forma, a lei da gravitação é 
consistente com o que observamos, mesmo com a isotropia do espaço violada, ou seja, as 
demais dimensões são distinguidas da dimensão extra pois as linhas de força só podem 
se propagar até uma distância re na dimensão compacta. Sendo a massa de Planck 
inversamente proporcional à Gy, temos em largas distâncias, com r > re, a adequação 
da lei da gravitação a um espaço-tempo 4-dimensional (3 dimensões espaciais e 1 di- 
mensão temporal) e sua diminuição com 1/(M,r)2. Em curtas distâncias, r < re, a lei 


n—2 


da força diminui com 1/(Mr) 


Até o momento, nos referimos a dimensões extras de tamanho muito reduzido, da 
ordem da escala de Planck. Para Teoria de Cordas, as 6 dimensões extras teriam esta 
ordem de grandeza. Entretanto, é possível pensar numa dimensão não tão reduzida 
e até infinita, como no modelo Randall-Sundrum. Restrições experimentais (veja [7]) 
indicam que tais dimensões extras podem ser menores que 10°GeV , isto é, 16 ordens de 
grandeza maior do que a escala de Planck, 10!ºGeV, e, além disso, através do trabalho 
de Arkani-Hamed, Dimopoulos e Dvali [8], que deu origem ao modelo ADD, a previsão 


para o raio da dimensão compacta no caso de 2 dimensões extras pode ser Imm. 


O modelo ADD utiliza um conceito proveniente da Teoria de Cordas que será modi- 


ficado e utilizado nos modelos Randall-Sundrum, o conceito de brana. O termo brana 
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deriva-se de membrana e significa uma superfície imersa num espaço com maior número 
de dimensões. Para Teoria de Cordas, existem subvariedades p+1 — dimensionais imersas 
no bulk 10-dimensional chamadas p—branas. Entre as p—branas existem as D—branas, 
que são extremidades de cordas abertas, 1—brana. Assim, no modelo ADD temos um 
espaço tempo 4 + d dimensional, com as d dimensões extras compactas, contornado 
por branas onde todos os campos, exceto o gravitacional, estão confinados nas subva- 
riedades. A motivação do modelo ADD foi resolver o problema da hierarquia, que é 
a diferença entre a escala de massa eletrofaca, Mp ~ 10°GeV, e a escala de massa de 
Planck 4-dimensional, M, ~ 10!°GeV. 


Entretanto, a dimensão extra pode ser maior do que alguns milímetros, pode ser 


infinita como veremos no modelo Randall-Sundrum 2 (RS-II). 


3.2 Modelos Randall-Sundrum 


Existem dois modelos Randall-Sundrum, ambos utilizam a idéia de um bulk de 5 di- 
mensões com todos campos, exceto o gravitacional, confinados na brana. O modelo RS-I 
consiste de duas branas 4-dimensionais imersas no bulk com o tamanho da dimensão 
extra finito. O modelo RS-II consiste de uma brana 4-dimensional imersa no bulk com 
o tamanho da dimensão extra podendo ser até infinito. Os modelos Randall-Sundrum 


foram desenvolvidos com o intuito de resolver o problema da hierarquia. 


Como veremos, devido ao fato de transportar energia, a gravidade será mais forte 
perto da brana e diminuirá exponencialmente quando nos afastarmos dela na dimensão 
extra. Dessa forma, a métrica na brana terá um fator exponencial que dependerá 
dessa nova dimensão e, assim, não será fatorável. Esse novo fator na métrica atenua a 


gravidade e corrige o problema da hierarquia. 


Randall-Sundrum I 


De acordo com [9], o modelo RS-I adota duas 3-branas imersas num bulk 5-dimensional 
anti-de Sitter, cuja coordenada extra, y, é compacta com raio re. Além disso, o bulk pos- 
sui simetria S! /Zo. Uma das branas pode ser identificada como nosso universo visível e 
a outra como um universo escondido. Ambas possuem posições fixas no bulk e podem 


ser interpretadas como condições de contorno neste espaço 5-dimensional. Assim, pode- 
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mos restringir o intervalo de interesse da coordenada y colocando uma brana em y = 0, 
a brana visível, e a outra em y = Tre, à brana escondida. Na figura 3.1, é representada 


uma imagem deste modelo com uma dimensão espacial suprimida. 


brana visivel 


PAN 5D brana escondida 


Figura 3.1: Representação do modelo RS-I onde uma dimensão espacial foi suprimida. As duas 
branas, visível e escondida, possuem tensão negativa e positiva, respectivamente. Ambas estão imersas 


num bulk anti-de Sitter com constante cosmológica Asp. 


Então, podemos escrever as métricas das duas 3-branas como 


Gois = G” (y = 0)ôpô,, (3.2) 
Jese = GY” (y = Tre) 515), (3.3) 


onde G”” é a métrica no bulk, os indices gregos vão de 0 a 4 e os índices latinos vão de 
Oa 3. 


A ação clássica que descreve este modelo é dada por 
S o Ssp + Svis + Sese (3.4) 


sendo S5D, Svis, Sese AS ações no bulk, brana visível e brana escondida, respectivamente. 


Podemos escrevê-las da seguinte forma 


Ssp = i J/|G|[-Asp + 2M, R]d zdy, (3.5) 
as = / V iza] eia Je 2rvis]0(y)d*xdy (3.6) 


Fo = É V [Jesc] (Lese P QresclO(Y Pos tr.)d'xdy, (3.7) 
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onde G, vis € Jese São OS determinantes das métricas do bulk, brana visível e brana 
escondida, respectivamente. Asp e Msp são a constante cosmológica e a massa de 
Planck em 5 dimensões, respectivamente. R é o escalar de Ricci, Avis é a tensão na 
brana visível e Aesc é a tensão na brana escondida, ambas atenuam o efeito da constante 
cosmológica do bulk na brana. As lagrangeanas L,;, € Lese das 3-branas não precisam 
ser conhecidas em detalhes para se determinar o métrica do bulk. As integrais acima 


são em relação a x”, coordenada na brana, e y, coordenada da dimensão extra. 


Podemos obter as equações de Einstein em 5 dimensões extremizando a ação (3.4). 


Dessa forma, as equações de campo são 


1 
WV _ yp ab cu 
y |G (r 5€ R) = a [2V/ |GvislGnisOs Op Avisd (Y) 
[GG Asp + 2y |Jesc|Jescðh 5p Aesc5(Y — Tre). (3.8) 


Assumindo que existe uma solução quadridimensional que obedece a invariância de 
Poincaré nas direções x”, quando a geometria é plana (Minkowski), pode-se usar a 
métrica 5-dimensional 

ds? = Wanda" da” + dy? (3.9) 


como Ansatz em (3.8). Assim, as equações de Einstein reduzem-se a 


A 
Go? = ——S5 (3.10) 
AMB, 
i waa Dede ) 
3 "”" vis yY | esc Y— Tle 11 
7 OME, 2MB, ely) 


(19 


Os símbolos linha indicam derivação com respeito a y. A solução para (3.10) que 


respeita a simetria Zo, y > —y, é 


Asp 


o = ly “AM? 


(3.12) 


Vemos em (3.12) que a solução aceitável é aquela com Asp negativa. Dessa forma, o 


espaço-tempo entre as duas branas é do tipo AdSs. 


Como a métrica (3.9) é uma função periódica em y, se derivarmos (3.12), válida no 
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intervalo —mr,. < Y < Tre, obtemos 


r l) Oly ate) fs (3.13) 


Então, obtemos (3.11) a partir de (3.13) se Avis, Aese € Asp estiverem relacionadas com 


um escalar k da seguinte maneira 
Avis = —Aese = FI2M5pk e Asp = —24MPk?. (3.14) 


As relações acima entre os termos cosmológicos do bulk e da borda são fundamentais 
para se obter uma solução que respeite a invariância de Poincaré 4-dimensional na 


brana. Com isso, a métrica do bulk é 


ds? = em dr dr? + dy). (3.15) 


Para ambos os modelos, RS-I e RS-II, a métrica do bulk é a mesma com a diferença 
de que para o modelo número um o sinal da exponencial é positivo e para o modelo 
número dois o sinal é negativo. Para o RS-I, o crescimento da exponencial na direção 


y é barrado pela brana escondida. 
A constante de Newton para a brana visível é dada por 


Nois 


Cine 
4D 6MS, , 


(3.16) 


onde a tensão negativa na brana visível induz uma gravidade repulsiva na mesma. Na 


brana escondida a tensão é positiva. 


Randall-Sundrum II 


No modelo Randall-Sundrum II, de acordo com [10], apenas uma 3-brana é utilizada 
e, como no RS-I, o bulk tem uma geometria anti-de Sitter, AdS;, com simetria Zə. A 
coordenada y pode tender ao infinito, ou seja, ao contrário do modelo anterior, o RS-II 


não utiliza uma dimensão extra compacta e, dessa forma, re — 00. 


A métrica do bulk no RS-II é a mesma do RS-I com sinal negativo na exponencial. 


A maneira que se obteve tal resultado é semelhante ao RS-I e, conseqiientemente, 
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omitiremos os cálculos. Assim, a métrica no bulk para esse modelo é 
ds? = e *lyda?da” + dy”. (3.17) 


A vantagem do RS-II reside no fato de que a única brana, a brana visível, tem tensão 


positiva, ou seja, possui uma gravidade atrativa. 
Fazendo k = 1/1 e introduzindo a coordenada z = le! em (3.17), obtemos 


[2 x 
ds? = (dz + nada" da”), (3.18) 
onde [é o raio AdS. A métrica de Minkowski em (3.18) pode ser substituída por qualquer 


métrica Ricci plana e será também solução de (3.8) como pode ser visto em [26]. 


3.3 Projeção das equações de Einstein na 3-brana 


A partir de agora, vamos obter as equações de Einstein na 3-brana (imersa num bulk 
5-dimensional com simetria Z2) com auxílio das equações de Gauss e Codacci, de acordo 
com [16]. Vamos, então, relacionar o tensor de Einstein na brana com elementos no 
bulk. Novamente, todos os campos, exceto o gravitacional, estão confinados na brana. 
As equações na 3-brana se reduzem às equações de Einstein convencionais no limite de 
baixas energias. Os índices gregos vão de 0 a 3 para tensores em 4 dimensões e de 0 a 
4 para tensores em 5 dimensões (a dimensionalidade será indicada pelos números entre 
parênteses inseridos no lado superior esquerdo do tensor). Esta seção foi diretamente 


influenciada pelo capítulo 6 de [27]. 


Nosso universo, uma brana 4-dimensional, será indicada por M e terá uma métrica 
dada por 


uv = Guv = Ny Ny e gee Ee qe as, nn”, (3.19) 


onde g, é a métrica no bulk e nº é um vetor unitário normal à brana. Podemos escrever 


algumas propriedades que nos ajudarão na obtenção das equações desejadas 


Vv V 


Q Ia TNN, mn =1 e qun” =0, (3.20) 


sendo q; um operador de projeção, isto é, dado um campo qualquer, X”, no bulk, q, 


calcula sua projeção na brana. 
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As equações de Gauss e Codacci podem ser escritas como 


CR ago O Roo lu Gs 0/05 + Ky" Kas — Ke" Ka (3.21) 


DK = DR Rana, (3.22) 


respectivamente, onde Kyy = q PV ong é a curvatura extrínseca de M, K,” = K 
é o traço de K,, e D, é o operador derivada covariante em relação à métrica quv- 


Contraindo (3.21) com relação aos índices a e y, obtemos 


O Rgs =® R" pol dg deeds” + K Kgs — K5° Ka, (3.23) 
como qe = g,° — Nn e qL = JP — NaN”, (3.23) fica 
Res = ent (g, = Nan Io o NaN” Jag 95° a K Kgs o K5” Kga 


= ORM olg — I Nan? — JPNN + nun Nan ag’ qs” + 


Be eek, (3.24) 
Usando as propriedades (3.20), (3.24) pode ser escrita como 


ORs = ORI Ger — Myung q + K Kos — K Kg, 
= Rega O RY poponun ggg + KK ps — Kj" Kga- 


Fazendo mudangas de indices, obtemos 


ORS?) Risada do —©) Re sro BN a a + KKm — DES Maga (3.25) 


Agora que temos o tensor de Ricci em 4 dimensões, podemos obter o escalar de 


Ricci contraindo (3.25) com a métrica q!” 
aR=— gh) Ruw, 
OR =©) Riog — gh” (5) RP nena, q +K? — K" Kpa- (3.26) 


ôpo 


Dessa forma, com auxílio de (3.25) e (3.26), as equações de Einstein ficam 
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1 
gim OR, 


Ried dy —©) R soon" 4,09," + KKw z Ky Kop = 


Og, = OR, — 


1 v o a 
= Mu Roo” = q" OR soca" q, Gy + K? =K“ Kual, 


após alguma manipulação algébrica, temos 


1 7 
er = MRO _ (5) Rood? 54w FG tg Ie Kua 


E ue — K” Kua — q!” ORP ong g q], (3.27) 


onde 


Ew =(5) RP soon Gy. (3.28) 


Como qu = sed a" o segundo termo do lado direito de (3.27) pode ser escrito da 


seguinte forma 


16) 1 


— Rie dw = odio O Riog”) 
= 5 (ate — nan”) ©) Rsolg™ — nn) qa” 
= [5 (06. nn CRY Ren nada, 
entao, 
Ra Oy = — Sir O Raa? + Randon qia, E 
d Sino eT le snanon'n® Rood a”. (3-29) 


Usando (3.20) e o fato de que OR = g% O Rso, (3.29) fica 


16) 


ES Ria quv = 


1 “ey SE a 
2 — 5960 © Rg? a, + 5 Poon'n quv- (3.30) 


2 


Assim, substituindo (3.30) em (3.27), obtemos 


1 1 ~ 
OG Es ( © Rso = 9950 ga + 9 luv © Renin? = Ew a 


Q 1 Q V o 
+ KK — K Kyo — 5d [K? — KK pa — gt” OR soon" Gu Gy], 
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sabendo que Gs = Rs, — 5950 (5) R, então 


DG a E GsoQ |? = Bey + FOR ap a DG Ties ate 


1 
+ zul Roonênº — K? + KK a + Re soon? q"). (3.31) 


Novamente, com auxílio de (3.19) podemos escrever (3.31) da seguinte forma 


DG = OG5oq,0q," — Ew + KKw — K a wot 


V 


1 
+ Sul Kpa — K? + Rion?n? + RY, nanl(g — nën”), 


2 
ou seja, 
OG = OGoogfa” — Ew + K Kw — Ky K uat 
1 
+ou |K" Ka — K? + Rona +) Rf ngn” —(5) Roron w n]: (3.32) 
Fazendo 
OR na =o" © Reign’ =" Ran? =E R nº (3.33) 
e 


ORF nanan? = 998 ©) Rgspongn n?n? =") Rgspon ninn = 0, (3.34) 


pois ©) Rgspo é um tensor anti-simétrico e nºnºn?nº é um tensor simétrico, o tensor 


4-dimensional em (3.32) fica 
DG = Gaga + E ze els t Kona EE 


1 
+Qu © Rsonên” — Iv SKK so): (3.35) 


Vamos, a partir de agora, relacionar MG uv Com O tensor energia-momento na brana. 


Para isso, vamos usar as equações de Einstein em 5 dimensões 
1 
Gio = © Rsc — a ORS Ta (3.36) 


onde «x2 é a constante gravitacional em 5 dimensões e “)75, é o tensor energia-momento 
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do bulk. Usaremos também o tensor de Weyl em 5 dimensões 
1 
BCs = “Reps — 3 Jor ©) Rsg — Gos ©) Rap — goy O Rio + 966 ©) Ra) + 


1 
a Tp (90456 — Jos 9) OR: (3.37) 


A partir de (3.36), obtemos, com auxílio de 9º”, 
(5) 2 2 (5) (5) KB fo (5) (5) 
R= ~ 3/5 T e Roo = 3 (3 Tso a. Tg) (3.38) 


ou seja, o escalar e o tensor de Ricci em 5 dimensões. Substituindo-os na equação 
(3.35), obtemos 


DG = ie Oia a = Ery + HE o Kee Ae 4 
k2 1 
f duz (8 Tio nn = 5d K? — K" Kpa). (3.39) 


Podemos relacionar Ery com Co usando (3.37). Sabemos que 


n b 6,0 
Ew = OR spo P BN Qp Wd > 
g” (5) Raspongn" q a,” 
= © Repo nado, 
depois de trocas de indices temos 
Biy) Rann nara, (3.40) 


Multiplicando (3.37) por n°n°q,°q,°, e com o auxílio de (3.40) e (3.20), teremos 


z 1 7 1 
Ew = 3! Regata’ +® Ryann Ie) = o © Rgspg, fa) + 


+ E Casta fa), 
e com ajuda de (3.38), É,» fica 


~ 1 K2 Ke a 
Ew = ala Tys —©) Tose)a fa) + ae OD — © Tyan 9659," 4,"] — 
1 2 

~ is © T9569,.°4,°) + Ep, 
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~ K2 1 1 a 
Ew = E Tag, q? = gh OT g¢s6q,°4,° a 3155 am nº quo + Ep, (3.41) 
onde 
Ew = seven ngs a. (3.42) 


Uma propriedade de E, é o fato de que seu traço é nulo, isto é, EH, = 0. Substituindo 
(3.41) em (3.39), obtemos as equações de Einstein 4-dimensionais em função do tensor 


energia-momento 5-dimensional 


X Ré he 
DG = =e? Tso hg? KK Ka = E Tosa qu? as E OUT Gis, JE 
4 Ks OT nn Qu + Ew) 4 5 (3 DT; —© Tgs,)n'n? — 
3 ya quv T uv) T luv 3 ôo Joo 
1 Q 
= PYAR KK (3.43) 


e, logo em seguida, fazendo uma troca apropriada de índices mudos, obtemos 


2 2 OT 
IG = E Tia. a + [9 Tsonên” ce -g lw) = Ew + 
1 
+ K Ky — Ki? Kia Qu — KY Kya): (3.44) 


2 


A expressão (3.44) nos revela que as equações de Einstein em 4 dimensões, nesse cenário, 
possuem termos que não dependem do tensor energia-momento. Dessa forma, tais 


termos fazem o papel de uma matéria geométrica que curva o espaço-tempo. 


A partir da equação (3.22), com o auxílio de (3.38) e do fato de que g;onº = 0, 


obtemos a relação 


2 
V K F 
D,K} — D,K = (3 OTe -® Tgso Pa, 


DKF = DK =n OT g (3.45) 


Até o presente momento, não utilizamos qualquer tipo de simetria em relação ao 
tensor energia-momento. A partir de agora, pensaremos num cenário brana e lançare- 
mos mão de algumas escolhas. Comecemos com a escolha de uma coordenada y, a 
coordenada da dimensão extra, e adotaremos uma única brana, que será localizada em 


y =0. Assim, n dx“ = dy e a métrica do bulk será dada por 


ds? = q,,dr"dzx” + dy”. (3.46) 
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Assumiremos que o tensor energia-momento no bulk tem a forma 

Ty = —Asp9Guv + Sw5(y), (3.47) 
onde 


Spo = Ada Tá (3.48) 


sendo Asp a constante cosmológica em 5 dimensões, À e Ty a tensão e o tensor energia- 


momento, respectivamente, na brana e Tn” = 0. 


O tensor energia-momento T, de acordo com (3.47) tem um comportamento singu- 
lar devido à brana, que divide o bulk em duas partes, quebrando a simetria de translação 
em y. Assim, existe uma região de transição entre o bulk e a brana e isso requer uma 


solução de continuidade que é dada pelas condições de junção de Israel [11]: 


E 1 
[Ky] = —K3 (Sv = 345), (3.50) 
onde o operador |X] corresponde a 
[Xv] = lim X- lim X = X*+- X7. (3.51) 
y—+0 vo 


A partir de agora, como nos modelos Randall-Sundrum, vamos impor a simetria 
Zə e fixar a brana num ponto fixo da dimensão extra. Com isso, a simetria imposta 


relaciona a curvatura extrínseca da brana com o seu tensor energia-momento 


£ 1 1 
Ko = -Kw EE 5 85 (Su z 34m9). (3.52) 
Por causa dessa simetria, podemos calcular coisas em ambos os lados, + e -, da brana. 


Substituindo (3.52) em (3.44) obtemos 


DG = —Napdu + 8TGNTu + Ket — Ey, (3.53) 
onde i i 
Aap = 5 5(Asp + gA”) (3.54) 
2) 
Gye (3.55) 


~ 48m 
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ge Etna? + ale + S guTagT” — ag ur”. (3.56) 

As equações (3.53) são, finalmente, as equações de Einstein projetadas na brana 

e se tornam as conhecidas equações da Relatividade Geral quando tomamos o limite 
de x — 0, mantendo Gy finito e desprezando E,,. O tensor E, é a parte do tensor 
5-dimensional de Weyl projetado na brana e carrega consigo a informação do campo 
gravitacional fora da brana e como ele a influencia. Portanto, E, depende da geometria 
do bulk e não pode ser determinado livremente pois é vinculado ao movimento da 
matéria na brana. Usando a identidade de Bianchi, D” OG u = 0, observamos a 


relação entre E, € Tw 


DYE, = K (DL Kop — DeKva), 
1 1 
= grat (Dr Tas — Dgtva) + 3 (Tuy — qu T) D”T]. (3.57) 


Vemos, então, que E, não é livremente determinado, mas sua divergência é vinculada 
à matéria na brana. Se decompusermos E,, em uma parte transversa, Es , € numa 
a a última será determinada completamente pela matéria. Dessa 


forma, se E: for zero, as equações (3.53) tratarão somente de elementos contidos 


parte longitudinal, E 


na brana pois a parte transversa corresponde a ondas gravitacionais no bulk. Assim, 
para resolvermos completamente as equações na brana deveremos, em geral, resolver as 


equações no bulk. 


No próximo capítulo, veremos como caracterizar e construir, de forma geral, soluções 


de buracos negros esfericamente simétricos na brana. 


Capitulo 4 


Buracos Negros Estaticos com 


Simetria Esférica 


A idéia de um objeto tao compacto e com um campo gravitacional suficientemente 
intenso de tal forma que a propria luz nao consiga escapar foi proposta originalmente 
por John Mitchel em 1783 [24]. O conceito de campo nao foi usado a principio, e toda 
argumentação foi embasada na Teoria da Gravitação Newtoniana. No século XX, a 
mesma idéia retornou com toda força depois da Teoria da Relatividade Geral e, devido 
a John Wheeler num artigo intitulado The Know and the Unknow, publicado em 1968 
na American Scientist e na American Scholar, um objeto deste tipo recebeu o nome de 


buraco negro. 


4.1 Horizontes 


Podemos representar um espaço-tempo de forma muito útil através dos diagramas 
de Penrose, pois regiões infinitas podem ser representadas por linhas finitas e pontos. 
Nestes diagramas, de forma pouco rigorosa, o infinito corresponde ao destino e/ou 
origem de curvas do tipo luz, tempo e espaço. Então, temos três tipos de infinitos: 
tipo luz, tipo tempo e tipo espaço, que correspondem ao destino e/ou origem de curvas 
do tipo luz, tempo e espaço, respectivamente (veja por exemplo |21|). Exemplo de um 
diagrama de Penrose para o espaço-tempo Schwarzschild anti-de Sitter, conforme [22], 


é mostrado na figura 4.1. 


36 
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singularidade 


singularidade 


Figura 4.1: O diagrama de Penrose do espaço-tempo Schwarzschild anti-de Sitter. No ponto A, 
podemos enviar sinais de luz que terminarao no infinito, r = oo, ou no horizonte de eventos. No ponto 


B, dentro do horizonte, estes sinais s6 poderao terminar na singularidade. 


No espaço-tempo de Minkowski, curvas do tipo luz terminam no infinito tipo luz. 
Entretanto, um espaço-tempo que contém um buraco negro possui uma região onde 
curvas do tipo luz não terminam no infinito tipo luz, curvas do tipo luz terminam numa 
singularidade, conforme é indicado na figura 4.1. Assim, o contorno desta região é 


chamado horizonte de eventos. 


Então, para definir um horizonte de eventos devemos conhecer a estrutura passada 
ou futura do infinito do tipo luz e, consequentemente, isso faz com que um horizonte 
de eventos seja um objeto não local. Ao contrário, um horizonte de Killing (veja [21]) é 
um objeto local, pois é definido como uma superfície H onde os vetores de Killing X” 
tangentes à ela são do tipo luz quando calculados em H, ou seja, X“X, |a= 0. Em [20] 
é mostrado que num espaço-tempo estacionário e assintoticamente plano o horizonte de 


eventos coincide com o horizonte de Killing. 


Métricas estáticas em 4 dimensões com simetria esférica podem ser escritas na forma 


geral 
2 


ds? = A(r)dt? — se r2(do? + sen?0d¢’). (4.1) 

Bir) 
A partir de um dado A(r), construiremos soluções estáticas esfericamente simétricas 
através do cálculo de uma família de funções B(r). Entretanto, antes de prosseguirmos, 
vamos formular as condições que a função A(r) conduz a métricas de buracos negros 


com horizontes de Killing. Para tal feito, utilizaremos a coordenada quase global, u, 
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que obedece a 


O motivo que nos faz adotar tal coordenada reside no fato de que (4.1) toma a forma da 
métrica de Schwarzschild e assim podemos utilizar toda ferramenta (extensão analítica 
da métrica, por exemplo) utilizada nessa solução bem conhecida. Utilizando a coorde- 
nada u, (4.1) fica 


du? 


ds? = A(u)dt? — AG 


— r?(u)(d6? + sen26dd), (4.3) 


onde 
Aw = A(r), raj=r-e Aw) (5) = B(r). (4.4) 


Nas geometrias em que A(r) = B(r) (por exemplo, solução de Schwarzschild), temos 


que u =r , ou seja, dr = du a partir de (4.4). 


Os zeros de A(u) representam pontos onde a métrica no sistema de coordenadas 
(t, u, 0, Q) é problemática. Uma condição necessária para que (4.3) descreva um espaço- 
tempo fisicamente razoável é que as funções A(u) e r(u) sejam CS (ou pelo menos 
diferenciáveis um certo número de vezes). Isto implica que em um zero, u = h, da 


função A(u) devemos ter 


Alu) ~ (u — h)* e r(u) ~r} +cte(u— h)º, (4.5) 


onde k é um inteiro positivo diferente de zero, s também é um inteiro positivo diferente 
de zero, que caracteriza o comportamento de r(u), r} é o raio do horizonte de eventos 
e o símbolo ~ relaciona grandezas de mesma magnitude e ordem dentro de um certo 
limite. Definimos a ordem do horizonte pelo número k. Quando k = 1, dizemos que o 


horizonte é simples; k = 2 denota um horizonte duplo e assim por diante. 


A seguir, vamos mostrar que a superfície dada por u = h é um horizonte de Killing. 
Como foi visto no capítulo 2, se uma métrica é independente de uma coordenada, o 
espaço-tempo descrito por ela tem um campo vetorial de Killing. A métrica (4.3) é 
independente da coordenada t, logo podemos escrever um campo vetorial de Killing 
como 


X= ] 
Ot’ (4.6) 
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ou expresso através de suas coordenadas como X” = (1,0,0,0), que satisfaz a condição 
LxGw =0. (4.7) 

A norma de um vetor desse campo é dada por 
X? = go X? X? = Au). (4.8) 


No horizonte de eventos, temos que A(h)=0 e esse vetor é do tipo luz pois sua norma 
é nula. Então, o campo vetorial dado por 0/Ot é um campo de Killing do tipo luz e 
a superfície em u = h é um horizonte de Killing, uma membrana de caminho único. 
Podemos, igualmente, mostrar que a métrica escrita na forma (4.1) tem um horizonte 
de Killing em r = r4, onde r, é raiz de A(r). Um campo vetorial de Killing pode ser 
escrito como (4.6), sendo suas coordenadas dadas por X# = (1,0,0,0), e sua norma 
dada por 

X? = go X? X? = A(r). (4.9) 


Em r = r4, temos vetores do tipo luz e, dessa forma, a superficie em r+} é um horizonte 
de Killing. 


Como pode ser visto por exemplo em [28], a definição de gravidade superficial, «+, 


em termos de um campo de Killing, X, é 
V(X#X,,) = 2X. (4.10) 
Tomando a derivada de Lie de k} em relação a X obtem-se 
Lx, =0, (4.11) 


ou seja, K+ é uma constante de movimento e, além disso, é constante no horizonte 
de eventos. Utilizando as coordenadas quase globais e um campo de Killing do tipo 
X” = (1,0,0,0), com auxílio de (4.8), (4.10) fica 


(4.12) 


A gravidade superficial, k}, tem uma analogia com um conceito da termodinâmica, a 
temperatura. Na termodinâmica de buracos negros, k, faz o papel de temperatura e, 
assim como na termodinâmica, a lei zero diz que a temeperatura, K+, é uniforme sobre 


um corpo em equilíbrio, ou seja, um buraco negro estacionário. 
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4.2 Extensão da métrica 


Veremos, a seguir, como obter a coordenada tartaruga, 7*, a partir da coordenada 
u, e estender analiticamente o espaço-tempo. A coordenada tartaruga em termos da 


coordenada r é = 
dr* 1 


Esso 4.1 
dr hfr)’ on 
sendo 
h(r) = y A(r)B(r). (4.14) 
Com auxílio de (4.4), podemos escrever 
A(r) 
dr = Do (4.15) 
e usar a relação acima em (4.13) para obter 
dr* 1 
a AD (4.16) 


Podemos utilizar (4.15), com auxílio de (4.14) e A(r) = A(u), na equação (4.12) para 
obter 


E 4.17 
k} == 
r=r4 
Definiremos, então, a coordenada avançada v 

v=t+f*, (4.18) 

onde neste novo sistema de coordenadas a métrica (4.3), utilizando (4.16), fica 
ds? = A(u)dv? — 2dv du — r (u)? (d8? — sen*6d¢’) , (4.19) 

e a coordenada retardada 

w=t-, (4.20) 

onde neste novo sistema de coordenadas a métrica (4.3), utilizando (4.16), fica 
ds? = A(u)dw* + 2dw du — r (u)? (do? — sen?6d¢") . (4.21) 


A métrica, escrita nas formas (4.19) ou (4.21), é completamente regular no horizonte, 
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onde A = 0. Dessa forma, a extensao da métrica para além do horizonte pode ser 
realizada, redefinindo os limites da coordenada u de forma conveniente. Em principio, 
o processo pode ser repetido mais vezes se necessário, até obtermos a extensão maximal 


do espaço-tempo. 


O comportamento da coordenada tartaruga, 7*(r), para valores grandes de r , como 
pode ser visto na figura 4.2, indica se a geometria estudada é assintoticamente plana, 
de Sitter ou anti-de Sitter, conforme (2.37), (2.38) e (2.39). Como dr*/dr é dada pelo 


inverso da função h(r), no caso de uma geometria assintoticamente plana 


1 
lim —~ =1 4.22 
rohr) (4.22) 
e a curva f*(r) tende a uma reta com inclinação de 45 graus em relação ao eixo r. No 
caso de uma geometria assintoticamente de Sitter, 7*(7) não está definida para valores 
de r maiores que um certo valor crítico, rc. Quando tomamos o limite para r > ro 


obtemos 


lim E (4.23) 
E, — E 
emo S 


e o gráfico de 7*(r) mostra que esta função diverge quando r = rc. Para o caso de uma 
geometria assintoticamente anti-de Sitter, 7*(7) tende a um valor constante quando r 
cresce, como consegiiência de 


1 
lim — = 0. (4.24) 


Nn 
j=) 


EN 
j=) 
7 


Ww 
oO 
T 


coordenada tartaruga 
N 
S 


pa 
j=) 
T 


Figura 4.2: O comportamento assintótico da coordenada tartaruga, 7*, é mostrado para as geometrias 
assintoticamente plana, A = 0, de Sitter, A > 0, e anti-de Sitter, A < 0. 
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4.3 Construção de soluções 


A partir da função A(r), podemos obter uma família de funções B(r) e cada mem- 
bro desta família corresponderá a uma métrica do tipo (4.1) diferente. Nesta seção, 


descreveremos o procedimento para obtenção dessas soluções. 


Podemos escrever os escalares de Ricci, R, e o de Kretschmann, «x, usando (4.1), da 


seguinte forma 


2 A A2 AB, 2/4 B 
e 1 = B = B rr Rx r | T r | T | É ba ` 
R al A 242 2AB r ( A B )| (825) 
e 
k = Ra" Roo!” =4Kj +8K; +8K3 +4K1, (4.26) 
onde 
g, — B (244 = A? A,B, 
4 A? AB) 
BA, B, 1-B 
ky = = Ki = ——. 
2 Or A’ 3 2r? É r2 


O símbolo subscrito r significa d/dr. A finitude de « é um critério de regularidade para 
geometrias baseadas em (4.1), pois este escalar é a soma das componentes do tensor de 


curvatura. 


Impondo que o lado direito de (4.25) seja uma função de r, R(r), e definindo 
f(r)=rB(r), (4.27) 
obtemos uma equação diferencial ordinária de primeira ordem 


r(2A(r)Am(r) — Ar) + 3A(r)A (r) py _ 2A(r)Q(r) 


f(r) ! A(r) g(r) f(r) i g(r) , (4.28) 
cuja solução geral é (veja [5]) 
= 2A(r)e” i r r eo r 
ro = | S seede], (4.29) 


onde 
g(r) = 4A(r)+rA,(r), (4.30) 
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= " A,(r) A 
r= f H (4.31) 


Q(r)=2-— r°R(r). (4.32) 


Dessa forma, escolhendo um A(r) liso e impondo um R(r), calculamos f(r) por (4.29) e 
obtemos imediatamente uma família de funções B(r) com o auxílio de (4.27). A métrica 
ficará completamente conhecida desde que a constante de integração C, proveniente de 


(4.29), seja fixada, ou seja, um membro da famíla de funções B(r) seja escolhido. 


Em universos brana, pode-se obter soluções das equações de Einstein projetadas 
na brana sem especificar o tensor E,,. Não discutiremos as propriedades do bulk, a 
existência de soluções neste espaço-tempo 5-dimensional é garantida por teoremas de 
imersão |23]. A partir das equações projetadas obtidas no capítulo 3 

D dos e v v 2V o mv 
G”, = Mpg”, +81GNT + KET”, — EY, (4.33) 
onde Asp, quv € Th são a constante cosmológica, a métrica e o tensor energia-momento 
na brana, respectivamente, Gy é a constante gravitacional, Ks é constante gravitacional 
em 5 dimensões e 7”, é dado por 
V 


Ra! 1 ag 1 2 
T u = qual Ü + TQ) + g dev TagT = gg T , (4.34) 


sabendo que T = 74, o traço de E”, é nulo, ou seja, E!, = 0 e que em 4 dimensões G”, 
é dado por 
V V 1 


podemos igualar (4.33) a (4.35), fazendo v = ju, e obter 
R = 4Ap — 8STGyNT + Kem. (4.36) 


Conseqiientemente, vemos que o escalar de Ricci R dependerá de variáveis definidas na 


brana, pois a constante gravitacional em 5 dimensões é dada por 


o 487ºGN 


T (4.37) 


2 
K5 


sendo A a tensão na brana. 


Em [5] é mostrado como a escolha do valor da constante de integração de (4.29) 


conduz a soluções diferentes de buracos negros estáticos e esfericamente simétricos na 
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brana. Para funções A(r) que podem ser escritas próximas ao horizonte na forma 
A(r) ~ (r — r,)!/8, onde s € N, e (4.30) seja maior do que zero na região de interesse, 
um buraco negro terá um horizonte simples se a constante de integração for maior do 
que zero e terá um horizonte de ordem 2 + ps, onde p é o expoente da aproximação de 
Q(r) = 2 — r2R(r) ~ (r — r,)? entorno do horizonte, caso a constante de integração 
seja nula. Para funções A(r) que podem ser escritas próximas ao horionte como A(r) ~ 
(r — r,)?/*, onde s € N, e (4.30) seja maior do que zero na região de interesse, um 
buraco negro terá um um horizonte de ordem 2 + ps caso a constante de integração seja 
nula e a função Q(r) se comporte como (r — r,)? próxima ao horizonte. Neste caso, 


está excluída uma constante de integração diferente de zero. 
Podemos encontrar as componentes do tensor energia-momento efetivo, T“”°F, na 
brana se escrevermos (3.53) como 


G”, + Nap = Gr, (4.38) 


sendo K? = 81G y = 87, de acordo com as unidades geométricas. Então, a densidade, 
ef ão radial, p?’ ão t ial, p‘!, efetivas sã t 
p% , pressão radial, pq: € a pressão tangencial, p , efetivas sao expressas em termos 


das funções A(r) e B(r) do seguinte modo 


PE BE 
= thin — eit = BS Be 
i de BH i, BA 

AND A RAD gai É Ap ? 


B |2A A A,B 2/A B 
0 o _ 2 ef _ rr ry f ry ry r 
GG STRE A Æ AB (i i E) (4:39) 


Numa brana onde há vácuo, 7!” = 0 mas T””®f pode não ser nulo, pois o tensor energia- 


momento efetivo contém a influência do bulk na brana. 


No próximo capítulo, utilizaremos o procedimento descrito anteriormente para obter 
métricas de buracos negros estáticos e esfericamente simétricos com constante cosmo- 


lógica diferente de zero na brana. 


Capitulo 5 
Soluções com constante cosmológica 


A introdução de uma constante cosmológica nas equações de Einstein muda drastica- 
mente o caráter assintótico das soluções esfericamente simétricas. Dependendo do sinal 
da constante cosmológica, temos uma geometria assintoticamente de Sitter ou anti-de 
Sitter, caso o sinal seja positivo ou negativo, respectivamente. Casos com a constante 


cosmológica nula já foram estudados por [5, 12]. 


5.1 Soluções assintoticamente de Sitter 
5.1.1 Raízes de A(r) e g(r) 


A partir de um A(r) do tipo Schwarzschild-de Sitter 
A(r) =1-—--—+ (5.1) 


onde m é a massa e a? = 12/R, sendo R o escalar de Ricci, utilizamos o procedimento 

descrito no capítulo anterior para a obtenção de uma família de funções B(r) na brana 
sem conteúdo de matéria. Podemos escrever (5.1) como 
ae 2 
SE aT — 2ma 

A(r) = (5.2) 


a2r 


e definir um polinômio para estudar as raízes de A(r). Assim, se escolhermos 
p(r) = =r? + a?r — 2ma’, (5.3) 
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então p(r), ou seja, A(r) terá um máximo em r = a/ v3 e um mínimo em r = —a/ v3. 


Queremos que A(r) tenha 3 raízes reais, e para que isso ocorra 


Da? 
p (<5) = sagt ame? >o (5.4) 


j (=) E E +2ma? <0. (5.5) 


Podemos ilustrar essa situação através do gráfico abaixo. 


p(r) 4 


=1 


V3 (3 


Figura 5.1: Gráfico de p(r) com duas raízes positivas e uma negativa. Os valores r = a/v3 e 


r= —a/v'3 indicam um máximo e um mínimo, respectivamente. 


Dessa forma, p(r) terá três raízes reais quando (5.4) e (5.5) forem satisfeitas si- 
multaneamente. O valor de a que satisfaz ambas condições é a > 3\/3m, e como 
p(0) = —2ma? < 0 e os valores de máximo e mínimo são eqiiidistantes da origem, te- 
remos duas raízes positivas, r4 e re, e uma raiz negativa, rn, como ilustrado no gráfico 
5.1. As raízes positivas têm um signficado físico; r4 é o raio do horizonte de eventos e 


re é o raio do horizonte cosmológico. Ambas são dadas por 


a a? + F(m,a)3 — iv3(a? — F(m,a)3) 
a VHF (ma) | ii 
i 1 |a? + F(m,a)3 + iv3(a? — F(m,a)3) 
eTa V3F(m,a)3 | a 
onde 


F(m,a) = a? (3v3m + V27m? = a?) ; (5.8) 
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Para obter uma família de fanções B(r), a partir do procedimento descrito no capí- 
tulo anterior, a função g(r) = 4A(r)+4A,(r) não deve ter zeros no intervalo de interesse 
estudado, r} < r < re, pois g(r) aparece no numerador em (4.29). A seguir, vamos 
mostrar que g(r) no caso Schwarzschild-de Sitter sempre tem para a > 3/3m uma raiz, 


rg, no intervalo de interesse. A função g(r) para esse caso é 


6m 6r? 
glr) =4- — -5 (5.9) 


r a 


e para estudar os seus zeros vamos reescrevê-la como 


= 6r? + 4a2r — 6ma? 


.1 
g(r) ar , (5.10) 
e determinar os zeros do polinômio 

g(r) = 6r? + 4a?r — 6ma?. (5.11) 


Para valores de a > 3\/3m, o polinômio q(r) terá três raízes reais, como pode ser visto 


a partir das condições 


2a 56a 
— |= —6ma? > 0 5.12 
2 (55) 273 o 
E 2 56a 
a a 
-— | = -—~ — 6ma’ < 0, 5.13 
( 55) 27/3 en 
sendo r = 2a/3/3 e r = —2a/3V3 pontos de máximo e mínimo, respectivamente, de 


q(r). As condições (5.12) e (5.13) serão satisfeitas para a > 81/3m/28. Como queremos 
três raízes reais para A(r), e para que isso ocorra a deve ser maior do que 3/3m, o valor 
de a que satisfaz as condições (5.12) e (5.13) não é incompatível com 3/3m. Assim, 
a > 3W3m é compatível com as condições que envolvem p(r) e g(r). O polinômio g(r) 
tem duas raízes para r positivo, pois q(0) = —6ma? e os pontos de máximo e mínimo 
são eqiiidistantes da origem. Logo em seguida, vamos mostrar que g(r) tem uma raiz 


entre as raízes positivas de A(r) definindo o polinômio 
d(r) = p(r) — ar), (5.14) 


que terá zeros onde A(r) = g(r). O polinômio d(r) tem dois zeros para r positivo quando 
a > 2v'5m, ou seja, não é incompatível com o valor adotado anteriormente, a > 3/3m. 
O fato de d(r) ter dois zeros para valores positivos de r, q(2a/3V/3) > p(a/v'3) e 


q(a//3) > 0 garante que q(r) tenha uma raiz entre as duas raízes positivas de p(r), ou 
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seja, g(r) tem uma raiz, rg, entre [r+, re]. Conseqiientemente, não existe um valor de a 
tal que g(r) tenha um zero fora do nosso intervalo de interesse e, devido a esse fato, o 


uso de (4.28) é problemático. A figura número 5.2 ilustra o que foi dito. 


— qr) 
— p(r) 


=1 


Figura 5.2: Se A(r) tiver 3 raízes reais, sendo duas positivas e uma negativa, a deve ser maior do que 
3v3m e, consequentemente, g(r), ou seja g(r), terá pelo menos uma raíz entre as duas raízes positivas 
de A(r) ou p(r). 


Entretanto, para a = 33m as funções A(r) e g(r) terão apenas uma raiz positiva e 
o valor desta raiz será o mesmo para ambas funções. Este resultado pode ser mostrado 


se escrevermos A(r) na forma 
A()=(r-ry(r—-ry)K(r), (5.15) 


onde K(r) não tem zeros em ry e re Com r} = re, A(r) e g(r) ficam 


D 
a 
3 
WN 
Il 


C-re PRP) 
A(r = ri K (r) + 7[2(r — ri) K(r) + (r —11)°K,(r)], (5.16) 


Q 
—— 
3 
sai 
| 


onde subscrito r indica a derivada em relação a r. Dessa forma, A(r) e g(r) terão 
zeros em r, e o intervalo de interesse agora é r} < r < oo. Assim, podemos utilizar 
o procedimento descrito no capítulo anterior e obter um resultado para (4.29) dentro 


desse intervalo. 
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5.1.2 Família de Soluções 


A função T(r), dada por (4.31), para este caso é 


go qn ees 


Maca 
sane a ie ei (5.17) 
3” 3m m|’ 


sendo r} = 3m com a escolha a = 3/3m. Obtivemos, em seguida, a função B(r) 
impondo R = 4A,p = 12/a?, de acordo com (4.36), com auxílio de (4.29) e (4.27). 


Com isso, a família de funções B(r) obtida é 


=1 


Aid ss aa 
v5 Hr +12) ¥8(r — r2)! + 


2r 2 


2m r? 
+[1— o DR (5.18) 


onde rı = 3m(v5 — 1)/2 e r2 = —3m(v5 + 1)/2 , somente definida no intervalo de 


interesse r} < r < œ quando a constante de integração, C, for nula. Para o caso que 


MAT m = 
oe _Cm(r + 6m) (=) a 


C # 0, a função B(r) não está definida dentro do intervalo escolhido e para C = 0, 
temos, de acordo com (5.1), A(r) = B(r), ou seja, a métrica Schwarzschild-de Sitter 


extrema (quando r, = re, veja [13]) já conhecida 
2 A 2 2m i E: Do Didi r2 eop 2042 
ds“ = bee ors dt — | 1 — — — —— |] dr*—r*d0°—r*sen*6d¢*. (5.19) 


Então, a familia de funções B(r) obtida só contém um membro, quando C = 0. 
5.2 Soluções assintoticamente anti-de Sitter 


5.2.1 Raízes de A(r) e g(r) 


A partir de um A(r) do tipo Schwarzschild anti-de Sitter 


A(r) =1-—+ — (5.20) 


CAPITULO 5. SOLUCOES COM CONSTANTE COSMOLOGICA 50 


onde, novamente, m é a massa e a? = —12/R, sendo R o escalar de Ricci, utilizamos 
o procedimento descrito no capítulo anterior para obtenção de uma família de funções 
B(r) na brana sem conteúdo de matéria. A função A(r) pode ser escrita como 

rê + a2r — 2ma? 


A(r) = (5.21) 


a2r 


Podemos estudar seus zeros, máximos e mínimos definindo o polinômio 
t(r) =r? + ar — 2ma?, (5.22) 
e calculando sua derivada em relação a 7 
t.(ry) = 3r? +a. (5.23) 


A equação (5.23) revela que (5.22) não tem máximos ou mínimos. Neste caso, A(r) 


tem 3 raízes: 1 raiz real positiva, que representa o raio do horizonte de eventos 
= (5.24) 


onde 


G(m,a) = à? (3/3m $ V27mê +a?) l (5.25) 


e duas raízes complexas conjugadas. No caso Schwarzschild anti-de Sitter, o nosso 


intervalo de interesse será r} < r < oo. 


A função g(r) = 4A(r) + rA,(r) neste caso fica 
g(r) = 4- — + — (5.26) 


ou 
_ 6r? + 4a?r — 6ma? 


a2r 


g(r) 


e podemos igualmente definir um polinômio para estudar suas raízes. O polinômio 


(5.27) 


v(r) = 6r? + 402r — 6ma? (5.28) 
e sua derivada em relação a r 


v,(r) = 18r? + 4a? (5.29) 
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serão usados para estudar os zeros da função g(r). Igualmente à função A(r), g(r) não 


tem máximos ou mínimos, tem uma raiz real positiva 


(5.30) 


Tg = 


onde 


H(m, a) = a? (32a? + 729m? — 27m) , (5.31) 


e duas raízes complexas conjugadas. Como g(r) = 3A(r) + 1 + 1?/a?, de acordo com 
(5.20) e (5.26), g(r) > A(r) e nos conduz ar, < r4}. Dessa forma, não teremos problema 
ao usar o procedimento descrito no capítulo anterior e obter um resultado para (4.29) 


dentro do intervalo r} < r < oo. À figura 5.3 ilustra o comportamento das funções 
A(r) e g(r). 


— A(r) 
| — g(r) 


aa 
~T 


Figura 5.3: A função g(r) não tem um zero no intervalo de interesse r+ < r < 00. 


5.2.2 Família de Soluções 


A função I(r), dada por (4.31), para este caso é 


2kı — cık 2 2 
rpo et O mi] DE 
4cy — ¢? 4cy — ¢? 4cy — c? 


— 1 2 
(r — cs) | pe ia | E tera | , (5.32) 
r4 — C3) 2 raro 


onde ky, k2, k3, ka, C1, C2 e c3 são constantes que dependem de m e a. Por simplici- 


E. in =] + taln | 
Py ( 
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dade, omitimos as expressões explícitas destas constantes em termos dos parâmetros da 
métrica. Impondo R = 444p = —12/a?, de acordo com (4.36), com auxílio de (4.29) e 
(4.27), a família de funções B(r) nesse caso só foi obtida numericamente. O gráfico de 


B(r) para vários valores da constante de integração de (4.29) é exibido logo abaixo. 


08 10 12 14 16 18 20 
rim 


Figura 5.4: Gráfico da função B(r) para vários valores de C com m = a = 1. O raio do horizonte 
de eventos e a raiz de g(r) são, com estas escolhas, r} =1 e rg = 0,8 , respectivamente. A curva 
azul representa o caso em que A(r) = B(r). As funções B(r), nos casos onde A(r) £ B(r), divergem 


quando r — rg, que é representada pela linha tracejada. 


O gráfico de B(r) mostra que para valores de C que tornam A(r) £ B(r), a função 
B(r) diverge quando r — rg. Isto pode ser mostrado se A(r), g(r), Q(r) = 2 — r2R, 


er) e eP() forem escritos em torno de r} como 
A(r) = A(rg) + Ol(r — ro)), (5.33) 
g(r) = G(rg)(r — rg) + Ol(r =r) (5.34) 
Q(r) = Q(r4) + Ol(r — rg) (5.35) 
eTO = TC) Orr] (5.36) 


eTO) = eM) + Ol(r — ro]. (5.37) 
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Substituindo as expansões acima na expressão para f(r), obtemos 


_ 2A(rje™ | f" Ne Olr\dr 
fe) = A | ate ota +c 


= Alr) eC — Q(ra)(r4 — 19] Er) + Ol — ro)" (5.38) 


Como eè") = 0, f(r) =rB(r) e 1/r =1/r, + Ol(r — r,)], chegamos à expressão para 


a função B(r) próxima a ry 


Sendo A(r,) < 0 e Q(r,) > 0, quando r — r, temos que B(r,) — œ e, conseqiiente- 
mente, para todos valores de C em que A(r) £ B(r), rg é o valor mínimo que coordenada 
r assume. Uma evidência de que o limite r — r, é de fato uma singularidade física 
é o comportamento do escalar de Kretschmann neste limite. Conforme pode ser visto 
na expressão (4.26), « diverge com r > rg . Ou seja, a singularidade é transladada de 


r — 0 para r > ry. 


Outro detalhe no gráfico de B(r), figura 5.4, é o fato que para r < r} as funções 
A(r) e B(r) têm sinais opostos. Este problema é uma conseqiiéncia do sistema de co- 
ordenadas adotado. Conforme visto na seção 4.2, as coordenadas (t,r,0, 4) são válidas 
somente com r > r}. Mudando-se para um sistema mais adequado, obtemos explicita- 


mente a métrica com a assinatura +2. 


O comportamento assintótico desta família de soluções foi observado através do 
gráfico A(r) — B(r) em função de r/m. Observa-se que para diferentes valores de C, 
sobretudo para valores que tornam A(r) £ B(r), a diferença entre estas duas funções 
tende a zero, ou seja, esta família de soluções tende à solução Schwarzschild anti-de 
Sitter, A(r) = B(r), quando r cresce, conforme a figura 5.5. Então, as soluções obtidas 
são assintoticamente anti-de Sitter, tendem à geometria (2.39). Este comportamento 
assintótico pode, também, ser observado com auxílio do gráfico da coordenada tartaruga 
em função de r. Conforme foi comentado no capítulo anterior, percebemos que esta 
família de soluções é assintoticamente anti-de Sitter independente do valor da constante 
C, ou seja, a coordenada tartaruga tende a um valor finito quando r cresce, de acordo 


com a figura 5.6. 
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5.x 104 


0 20 40 60 80 100 


Figura 5.5: Gráfico da difrença A(r) — B(r) em função de r/m para valores de C que tornam 
A(r) Æ B(r). Neste gráfico foi usado m = a = 1. 
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Figura 5.6: Gráfico da coordenada tartaruga, 7*, em função de r. A curva azul representa o caso 
em que A(r) = B(r). Adotamos para construir este gráfico m = a = 1 e, conseqiientemente, r+ = 1, 
indicado pela linha pontilhada. 


Podemos obter o tensor energia-momento efetivo na brana para diferentes valores 


de C com auxílio de (4.39). Com isso, a densidade efetiva, pressões radial e tangencial 
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podem ser escritas como 
pf G‘, + 4A4p 
= DD, 


= (5.39) 
G", + 4A 
ef _ r 4D 
=" A 
Prad 87 (5 0) 
° G? A 
ef Gg + 4Aap 
= —— 5.41 
Pi Em ( ) 


O comportamento das componentes do tensor energia-momento efetivo pode ser avali- 
ado através dos gráficos abaixo e nota-se que conforme r cresce as componentes desse 
tensor tendem a um valor fixo, independente da constante C. Tal valor fixo corresponde 
ao caso onde A(r) = B(r), ou seja, a métrica Schwarzschild anti-de Sitter. Embora este- 
jamos tratando de uma brana sem matéria, temos dois ingredientes a mais: a constante 
cosmológia e o bulk. Dessa forma, as componentes do tensor energia-momento efetivo 


não são nulas. 


p efetiva 


08.10 12 14 42416 42418 20 
rim 


Figura 5.7: Gráfico da densidade efetiva, pº!, em função de r/m. A curva azul é o caso onde 
A(r) = B(r). Adotamos para construção deste gráfico m = a = 1 e o raio do horizonte de eventos, 
consegiientemente, é r} = 1. A densidade efetiva diverge, para os casos onde A(r) Æ B(r), quando 


r>rg= 08. 
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Figura 5.8: Gráfico da pressão radial efetiva, p®!,, em função de r/m. A curva azul é o caso onde 
A(r) = B(r). Adotamos para construção deste gráfico m = a = 1 e o raio do horizonte de eventos, 
conseqiientemente, é r, = 1. A pressão radial efetiva diverge, para os casos onde A(r) Æ B(r), quando 


r— ry = 08. 
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Figura 5.9: Gráfico da pressão tangencial efetiva, pý, em função de r/m. A curva azul é o caso onde 
A(r) = B(r). Adotamos para construção deste gráfico m = a = 1 e o raio do horizonte de eventos, 
conseqiientemente, é r+ = 1. A pressão tangencial efetiva diverge, para os casos onde A(r) 4 B(r), 


quando r — rg = 0.8. 
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5.2.3 Comportamento próximo ao horizonte 


Agora estamos interessados no comportamento desta família de soluções próximo 


ao horizonte de eventos, r+. Conseguiremos isto escrevendo (5.20) como 
A(r) = Ary)(r — r4) + Oltr — 1)?). (5.42) 
Além disso, A, e g(r) podem ser escritas como 


A, = A,(r4) Ar Apr (r)(r = r+) ae Ol(r a r+) (5.43) 


4A(r) + rA,(r) 
= Ari), + Ol(r —r4)]. (5.44) 


g(r) 


A função I, dada por (4.31), foi calculada próxima ao horizonte, e as exponenciais, 


usadas no cálculo de B(r), foram escritas como 


eT —14+O0[(r —1,)] (5.45) 
e = 1+ Ofļ(r — r4). (5.46) 

A função Q(r) fica 
Q(r) =2 -= r°R = Q(r4) + Ol(r—ry)], (5-47) 


e 1/g(r)? pode ser escrita como 


1 1 
g(r)? z A,(r4.)?r2 Ol(r — r+)]- (5.48) 


Conseqiientemente, B(r), dada por (4.27) e (4.29), pôde ser escrita como 


2Ae*! [ |” E 
Bj = rae? if g(r)Q(r)je “dr +C 
2C 2Q(r+) 


= Gan” = Rd A (r = r4) + Ol(r — r4)°], (5.49) 


próxima do horizonte. Dessa forma, uma família de funções B(r) foi calculada ana- 
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liticamente em torno de r,, e para a constante de integração, C, diferente de zero 


temos 2C 
= = ces RO 
B(r) = ECA (r— r,)+Ol(r — r4)'], (5.50) 
onde para C = Co, com 
1 
Co = zA ri (5.51) 


temos A(r) = B(r), ou seja, a métrica Schwarzschild anti-de Sitter. O resultado em 
(5.51) foi testado extensivamente com as ferramentas numéricas desenvolvidas. Nos 
gráficos apresentados nesta dissertação, fazendo m = a = 1 e, por conseqüência, r} = 1, 
o valor de C que torna A(r) = B(r) é C = Co = 8. Para o caso em que Č = 0, a função 
B(r) próxima ao horizonte será 

2Q(r+) 


B(r) = r? CHa Olea"). (5.52) 


Podemos escrever a função h(r), dada por (4.14), quando C ¥ 0, com auxílio de 
auxilio (5.42) e (5.50), como 


hr) = -rOle rs) (5.53) 


e a gravidade superficial, dada por (4.17), neste caso fica 


| C 


Para o caso em que Č = 0, conhecido como caso extremo, usamos (5.42) e (5.52) para 


obter 


hlr) = PED, rg)? + Olr = rt], (5.55) 


T+ 


sendo a gravidade superficial nula neste caso 
Em geral, com C # 0, a gravidade superficial é uma função de m,a e C, dada por 


1/2 


(5.57) 


See | 3V/3G(m, a)C 
ii 2(G(m, a} — a2)? 


CAPITULO 5. SOLUCOES COM CONSTANTE COSMOLOGICA 59 


com G(m,a) apresentada em (5.25). 


Como foi dito no capítulo anterior, existe uma analogia entre a gravidade superficial 
e a temperatura na termodinâmica. A lei zero da termodinâmica estabelece que a tem- 
peratura de um corpo em equilíbrio é constante, da mesma forma que kK, é constante 
no horizonte de eventos de um buraco negro estacionário. Buracos negros extremos são 
interessantes do ponto de vista da terceira lei generalizada da termodinâmica. Esta lei 
afirma que buracos negros com K+ = 0 não poderiam ser gerados por processos dinâmi- 
cos a partir de buracos negros com K+ + 0 [28]. O gráfico da gravidade superficial, para 


determinados valores de C, em função de m e a é exibido na figura 5.10. 


Figura 5.10: Gráfico da gravidade superficial em função da massa m e a. Utilizamos C = 2. 


Através de (4.4), podemos obter a coordenada quase global u(r) e calcular a função 
A(u), próxima ar,, para determinar a ordem do horizonte de eventos. Para o caso em 


que C > 0, temos um r(u) na forma 


1 2C 
rere + A, (rs) a = h), (5.58) 
que nos conduz a 
2C 
Alu) = [Su h), (5.59) 
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com um horizonte simples e com a constante de proporcionalidade igual a 2k}. Para o 


caso C = 0, temos 
Q(rs)? 


maja (5.60) 


r xr + 


com A(u) indicando que o horizonte de eventos, u = h, neste caso tem ordem 2 


Rij a eye (5.61) 


2A,(r4)r4 


Capitulo 6 
Perturbacoes Escalares 


Neste capítulo estudamos a estabilidade à evolução de um campo escalar da familia 
de soluções obtida no capítulo anterior. Os pioneiros no estudo sobre estabilidade de 
buracos negros foram Regge e Wheeler quando em 1957 no artigo [30] estudaram a 
solução de Schwarzschild e mostraram sua estabilidade. É interessante considerar a 
dinâmica de campos para o caso de espaços-tempos assintoticamente anti-de Sitter. 
Além de aspectos ligados à evolução do campo teste, a descoberta da correspondência 


anti-de Sitter /teoria de campo conforme (AdS/CFT) torna a investigação mais atraente. 


6.1 Introduzindo o campo escalar 


Em geometrias esfericamente simétricas 


dr? 
Btr) 


ds? = A(r)dt? — — r*(do + sen”6d¢’), (6.1) 
a propagação de campos escalares pode ser tratada com bastante generalidade, se con- 
siderarmos o background como fixo, e, portanto, o tensor energia-momento associado 
ao campo escalar desprezível. À equação de Klein-Gordon para campos escalares não 


massivos, ®, é dada por 


&=0, (6.2) 


onde [O é o operador d' Alembertiano 


= 9” VV, (6.3) 
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A equação (6.3) pode ser escrita em termos das coordenadas esféricas (t, r, 0, œ) 


= =) (V=39"9,), (6.4) 


onde g = det (g u). Como a matriz (g,,) é diagonal, o seu determinante é simplesmente 


o produto dos seus elementos da diagonal principal, e, portanto, 


V-9=4| Do) r2send. (6.5) 


Podemos escrever (6.4) em duas partes 


= Lh + Vo, (6.6) 


onde os operadores Ll e Və são escritos como 


a = [o (V=9906) + à (/=99"21)] (6.7) 


Vo = = [ð> (W/—9 97702) + Os (Vg g”03)] . (6.8) 


Substituindo os elementos da métrica (6.1) em Lz e V2, obtemos, após os devidos 


cancelamentos, 
1 8 8? 2 | A,(r) | B,(r) 9 
2 A(r) Ot? 7 aUo E | 2A(r) ! a Bl) ar (6.9) 
onde subscrito r indica d/dr, e 
1 1 0 ð 1 2 
V2 = ES [= (seno) O al . (6.10) 


Vemos, então, que o operador V2 é o operador Laplaciano em 2 dimensões escrito em 


coordenadas esféricas. 


Como o operador Ll só atua nas variáveis t e r, enquanto que o operador V2 só 


atua nas variáveis angulares, podemos expandir o campo & em multipolos, ou seja, 


p = X delt, r)Yem(9, ¢), (6.11) 
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onde Y,m sao os harmônicos esféricos, auto-funções do operador Laplaciano bidimensi- 
onal i 
V2Y em = matt +1)Yem, (6.12) 


soluções de 


LL +1)Ye m = 0. (6.13) 


1 ð A gorem | | Yem 
send 06 99 “sento 04? 


Substituindo o Ansatz (6.11) e (6.12) na equação de movimento (6.2), obtemos 


XO {0 (We) Yem + be Va Y em) = 0 
tm 


t(t+1 
eee À 2 (we) + 72 eh = 0. (6.14) 
tm 
Como as funções Y,m formam um conjunto completo, cada termo da somatória anterior 


deve se anular. Assim, ficamos com equações independentes para cada modo £ 


2 (We) + Dy, =0. (6.15) 


Com auxílio de (6.9), a equação acima fica 


Pe ap (Zee) E É = ay ol ESSA 


(6.16) 


ot? Or? r 24 2B Or a r2 


Temos então para cada modo ¢ uma equação parcial linear hiperbólica para a função 


pelt, r), que descreve a dinâmica da componente / do campo escalar ®. 


Para resolvermos a equação (6.16) numericamente, substituimos a coordenada r pela 


coordenada tartaruga 7* e definimos 
W(r,t) = 9(r,t)e(r), (6.17) 
onde o subscrito / será abandonado a partir de agora, 


a(r) = A(r)B(r) (6.18) 


RES: (6.19) 
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Depois de substituir as relações acima em (6.16), obtemos 


Oy Or O 2 10704 1 OF OB 
Pa — “ar Bp (PP) E "24 Or OF" 2B Or | 
or o |or* O 

Cor or E or 


(20) = —A [CM] gp. (6.20) 


Parte do segundo termo do lado esquerdo da equação acima pode ser escrito como 


2 1 Or* OA | 1 Or* OB _ 2 Lor 104 10B 
r 2A ðr OF* 2B Or OF* r 20r À4Ad Bom 
2 LOPE Yao 
~ r 20 (Zna | mB) 
2 1 o ða 
~ r Qa Or OF GE 
Substituindo (6.21) na equação (6.20) dividida por p, temos 
Pp adr Op, 09 \ (2.5. Lorde), - 
o? por Por lar) |r 2a ar om 
_ aor" Op | Op\ O (OR | 
p Or (“om | Pam) am \ Or 
a (OFN? / &p 79% 00, PN e(€+1) 
ce (es ae Bs 05) E -a| r2 le O23) 
Podemos escrever a equação acima reagrupando termos da seguinte forma 
Old o. OR Op si) SOR Or da OM OF CORY. 
o? Var) om (ar lr'zadror)' Or OF or) 
20 [Or * Ap Op adr Op [2 1 OF da) adr ð (ƏN Op | 
p \ðrj a| OF p Or OF kr ` 2a Or OF) p Or OF* \ Or j or’ 
a (or 2 92p LL +1) 
(5) El o= -ae = |. (6.23) 
Para que a equação acima fique do tipo 
O(r,t Op(r,t 
ye, (6.24) 


YOR? 


Ot? 


devemos impor 


ù (=) aif (6.25) 
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ay Pr sa EA o 
(ELE E (E) (E Bco q 
Vert =a |] -vo (627) 
onde 
ok ~x ~x ~x ae \ 2 

V0) =o eae (a tana an) te orar (or) arta (ar) ae (628 

Portanto, a equação (6.25) fica 

dr* 1 1 


= (6.29) 


que recupera a definição apresentada em (4.13). A partir de agora usaremos a derivada 


d e não mais a derivada parcial O. Logo em seguida, substituindo (6.29) em (6.26), 


obtemos 
2 lda l da 2dp | 
r 2adr 2adr pdr 
dp p 1 
dr T P r ( ) 


Assim, a função de onda (6.11), que descreve o campo escalar não massivo, fica 


pelr 


®(r,t,0,¢) = Dy, mb 6). (6.31) 


Substituindo a expressão para p na equação (6.28) e com auxílio de (6.29), chegamos à 


forma final do potencial 


1 da 
VS eae (6.32) 
e do potencial efetivo 
LL 1) 1 da 

Ver) =A | 2 | da (6.33) 

que pode ser escrito como 

t(L+1 1 

Ves(r) = A | ( 72 | | E [A(r)B,(r) + B(r)A,(r)] . (6.34) 


Pelo que foi visto no capítulo anterior, V., é positivo em ry < r < ov, pois A(r) 
e B(r) são funções positivas monotonicamente crescentes neste intervalo (veja figura 


5.4). Este é um resultado importante, pois a positividade do potencial efetivo garante 
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a estabilidade do campo escalar, conforme pode ser visto em [31]. 


Vamos estudar o comportamento de V.;(7*) próximo ao horizonte para a familia de 
soluções assintoticamente anti-de Sitter obtida no capítulo anterior. Para o caso em 


que C # 0, a função h(r) em primeira ordem de aproximação, de acordo com (5.53), é 


dada por 
hlr) = 2ny(r— r4) + O[(r — r4)?] (6.35) 

e a coordenada tartaruga, 7*(r), definida em (4.13), torna-se 

~ 1 1 

P(r) = —In(r > r4) + Of|(r — r,)2], (6.36) 

264 

ou seja, 

(P) = r4 +e (1+0(r-r5)3]). (6.37) 


Então, o potencial dado por (6.34), com auxílio de (5.42) e (5.50), quando £ > 0, será 
Ves(r)  (r=— 14) + Ol(r — ry", (6.38) 

e com auxílio de (6.37) fica 
Ves (#*) = Veg(r(#*)) ~ eletro, (6.39) 


Quando nos aproximamos do horizonte, 7* — —oo e o potencial decai exponencialmente. 


Para o caso em que C = 0, a função h(r), de acordo com (5.55), pode ser escrita 


como 


_ AGA) 


T+ 


h(r) (=r)? +Ol(r —14)?], (6.40) 


e a coordenada tartaruga 7*(r) torna-se 


2 1 
P(r) = — dis (r=r,)- 3. (6.41) 
V2Q(r+)Ar(r4) 
Como 7* < 0, próximo ao horizonte, podemos escrever (6.41) na seguinte forma 
2 1 
(r-r)? =- i =, (6.42) 
V2Q(r+)A (r+) f 


ou seja, 


(6.43) 
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Assim como no caso C + 0, utilizando (5.42) e (5.52), o potencial pode ser escrito 
como (6.38) quando £ > 0, e, com auxílio de (6.43), pode ser expresso como função da 
coordenada tartaruga 

Ves (7") ~ se (6.44) 


pre 


Quando ¢ = 0, 0 potencial efetivo sera 
Ves(r) ~ tr — ra)? + Ol(r —r4)°], (6.45) 
e pode, igualmente, ser escrito como função de 7* 
e 1 
Ve) ae (6.46) 


Então, quando Č = 0, o potencial efetivo tende a zero como uma lei de potência próximo 


ao horizonte de eventos. 


Podemos escrever (6.31) de tal forma que a parte temporal seja expressa através de 


uma exponencial, ou seja, 


&(r,t,0,6) = >> any, meee”, (6.47) 
Lym 


r 


e, conseqüentemente, (6.24) toma a forma de uma equação do tipo Schröedinger 


saz + (w? — Veslr)) e(r) =0, (6.48) 


que é a equação que governa a evolução do campo escalar não massivo para um dado 


modo £. 


O problema que estamos tratando, a evolução do campo escalar não massivo, é 
parecido com o problema de ondas eletromagnéticas emitidas por uma fonte numa 
cavidade, sendo este fenômeno governado por uma equação do tipo (6.24). A geometria 
da cavidade impõe uma condição de contorno ao problema fazendo com que ondas 
eletromagnéticas existam somente dentro. Para uma geometria assintoticamente plana, 
como Schwarzschild, o potencial cai como uma lei de potência com r — oo e existem 
soluções de (6.48) para cada w que correspondem a uma onda originária do infinito, 
sendo parte absorvida e parte refletida pelo buraco negro. As freqüências de oscilação, 
w, são chamadas freqiiéncias quase normais e são determinadas pelo potencial Vef, que é 


determinado pela geometria do espaço-tempo. Conseqüentemente, as freqüências quase 
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normais podem indicar a estrutura do espaço-tempo, inclusive a existência de dimensões 
extras como propõe Maartens et al. em [29], da mesma forma que as frequências das 
ondas eletromagnéticas no problema da cavidade são determinadas por parâmetros da 
própria cavidade. Numa geometria assintoticamente plana, as condições de contorno 
impostas são 

pm e” quando 7* — —oo (6.49) 


yet quando 7* — oo. (6.50) 


Para geometrias assintoticamente anti-de Sitter, como a família de soluções obtida no 
capítulo anterior, o potencial efetivo diverge no infinito e uma condição de contorno 


que pode ser imposta é 
w=0 quando f* =0, (6.51) 


como pode ser visto em [32, 31]. 


6.2 Cálculo numérico 


Será usado o chamado “problema das duplas coordenadas nulas” [33]. Na versão 
(1+1)-dimensional deste esquema, a equação hiperbólica de movimento é escrita usando 


as coordenadas nulas w e v, definidas no capítulo 4. Através destas coordenadas, vemos 


que 
o? o? o? o? 
o ðw? Ou. OU (P2 
e 
o? o? o? o? 
E e a A 
Se substituirmos as relações acima na equação de movimento (6.24), teremos 
ylw, v) 1 
Bud E -3V (rw, vlw, v) ; (6.54) 


O campo, então, é especificado na fronteira de um ângulo delimitado pelas semi-retas 
w = wo(v > vo) e v = vo (w > wo), que se interceptam no ponto (wo, vo), conforme 


ilustrado na figura 6.1. 
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Figura 6.1: Ilustração do problema das coordenadas duplas. 


Na maior parte das simulações feitas, usamos como perturbação inicial um pulso 


gaussiano no eixo w = wo, e fizemos o campo nulo no eixo v = Up, isto é, 
p(w, vo) =0, (6.55) 


wadi, 


bi) = exp | US 


(6.56) 


onde o é um parâmetro da gaussiana. Testamos diretamente se existe alguma de- 
pendência na forma particular do pulso através do uso de outras configurações iniciais 
alternativas. A menos de particularidades na região transiente, os resultados não apre- 


sentaram mudanças significativas. 


A partir do campo especificado na fronteira, o campo no interior do ângulo é deter- 
minado de maneira única. Então, é obtida uma forma discretizada para a equação de 


movimento (6.54), sugerida por [33], escrita como 


V(C)AwAv 


= (WB) + WD) +O(AS, (657) 


WA) = Y(B) + Y(D) — ¥(C) — 


sendo 
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= v(B)—-v(C) = O(A). (6.58) 


De maneira resumida, o algoritmo numérico consiste na varredura do plano w — v 
discretizado, onde dado três pontos calcula-se um quarto, usando a expressão (6.57). 


A integração começa em (vo, Wo) e prossegue no sentido de w crescente. 


6.3 Resultados numéricos 


O resultado obtido analiticamente deixa claro que o potencial efetivo em função da 
coordenada tartaruga, V.y(7*), próximo ao horizonte, tende a zero exponencialmente ou 
como lei de potência, de acordo com o valor da constante C. Este mesmo resultado foi 
observado diretamente através do cálculo numérico, como pode ser ilustrado através do 
gráfico semi-log! 6.2 onde temos um decaimento exponencial do potencial efetivo para 


valores de C maiores do que zero. 


Figura 6.2: Gráfico semi-log do potencial efetivo próximo do horizonte de eventos, para vários valores 
de C com C > 0. As linhas retas neste gráfico indicam o decaimento tipo exponencial. Os parâmetros 
do gráfico são m=1,a=lel=O0. 


'Logaritmo de base 10. 
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Para o caso extremo, C = 0, o decaimento do potencial efetivo é do tipo lei de 
potência, como pode ser visto no gráfico 6.3, sendo mais evidente quando ¢ = 0, pois, 
como vimos anteriormente, para este caso o potencial decai com 1/7** e para £ £ 0 


decai com 1/72. 
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Figura 6.3: Gráfico log-log do potencial efetivo próximo do horizonte de eventos no caso extremo, 
para vários valores de /. As linhas retas neste gráfico indicam o decaimento tipo lei-de-potência. É 
evidenciada a diferença entre a curva com £ = 0 e as curvas com £ > 0. Os parâmetros do gráfico são 
m=l,a=leC=0. 


Observamos, pelos resultados numéricos obtidos, que o campo escalar não massivo 


tende a zero quando t cresce e é limitado, ou seja, 


p(t, | < Ymar- (6.59) 


Além disso, o tipo de decaimento dependerá de C, a constante de integração de (4.29). 
O decaimento do campo possui duas fases: 

e Fase transiente: fortemente dependente da perturbação inicial; 

e Fase assintótica: função apenas dos parâmetros da geometria (m, a e C). 
Estaremos interessados na região assintótica. A figura abaixo ilustra o compartamento 


do campo em função do tempo, a transição entre essas duas fases é indicada pela linha 


tracejada. 
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Figura 6.4: Fases transientes e assintóticas do decaimento do campo escalar não-massivo, para um 
valor fixo de r. As figuras superior e inferior mostram o mesmo gráfico, em escalas normal e semi-log, 


respectivamente. A linha pontilhada separa a fase transiente da fase assintótica. Os parâmetros do 


gráfico são m = 1, a = 1, C = 10 e £= 1. 


6.3.1 Decaimento oscilatório-exponencial 


Para valores de C maiores que um certo valor crítico, Ceri, O campo é dominado 


pelo decaimento oscilatório amortecido exponencialmente 
w By emt = Er (6.60) 
onde w é um número complexo dado por 
W = tr + iwi, (6.61) 
sendo w, e w; a parte real e imaginária de w, respectivamente. Verifica-se que 


w; <0, (6.62) 
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e, dessa forma, o campo decai a zero exponencialmente. Com a massa e a constante 
cosmológica fixas, a medida que C diminui (com C > Ceri), O campo decai mais 
rapidamente, ou seja, |w;| aumenta. Esse ponto pode ser melhor ilustrado na figura a 


seguir. 


C=100 


ern 
= 
Sa 


m 
> 


Figura 6.5: Gráficos semi-log ilustram o decaimento do campo escalar em sua fase oscilatória- 
exponencial para dois valores diferentes de C. Os parâmetros do gráfico são m=1,a=1lel=0. 


6.3.2 Decaimento exponencial não-oscilatório 


Para C < Cerit, O comportamento do campo muda e o seu decaimento é exponencial 
sem oscilações, ou seja, 


Y ~ e (6.63) 


com k < 0. Com a massa e a constante cosmológica fixas, a medida que C aumenta 
(com C < Cerit), O campo decai mais rapidamente, ou seja, |k| aumenta. Ilustramos 


esse ponto na figura a seguir para dois valores distintos de C. 
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Figura 6.6: Gráficos semi-log ilustram o decaimento do campo escalar em sua fase não-oscilatória- 
exponencial para dois valores diferentes de C. Os parâmetros do gráfico são m=1,a=1lel=0. 


Os resultados numéricos sugerem que a transição do decaimento não-oscilatório para 
o decaimento oscilatório ocorre quando |k| > |w;|. Este comportamento foi observado 


para todos os valores de / considerados e essa transição é ilustrada na figura abaixo. 
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Figura 6.7: Transição entre o decaimento não-oscilatório para o oscilatório. A linha pontilhada 
marca o valor aproximado de Cerit, |k| = |w;| = |Imw|. Os parâmetros do gráfico são m = 1, a = 1 e 
L=0. 
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Podemos ter uma idéia do valor aproximado de Ceri com auxílio de gráficos semi-log 
do decaimento do campo escalar para alguns valores de C próximos de Cerit e observar 


a mudança no tipo de decaimento. 
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Figura 6.8: Gráficos semi-log do campo escalar próximo do valor de C = Cerit = 5,2 onde o 
decaimento passa de exponencial não-oscilatório para oscilatório. Os parâmetros do gráfico são m = 1, 
a=lel=0. 


Por exemplo, param = a = 1 e £ = 0, pode-ser notar que em Certi © 5,2 existe uma 
mudança no tipo de decaimento: para Certi < 5,2, o campo escalar não massivo decai 
exponecialmente sem oscilar e para Ceri > 5,2 o campo decai exponencialmente com 
oscilação. O fato de não termos valores exatos de Cerit é devido à limitação do método 


numérico. 


6.3.3 Decaimento em lei de potência 


Para C = 0, temos novamente uma mudança qualitativa no decaimento do campo 
escalar. Vemos que neste caso o decaimento do campo é como lei de potência. A 
mudança é esperada, visto que o comportamento do potencial efetivo com 7*, no limite 


de 7* — —oo, também muda com C = 0. Ilustramos esse ponto na figura 6.9. 
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Figura 6.9: Gráfico log-log mostra o decaimento do campo escalar não-massivo no caso extremo 
(C = 0). Neste gráfico, a linha reta indica o decaimento tipo lei-de-poténcia. Os parâmetros do gráfico 
são m = 1,a=1e£4=1. 


Um decaimento deste tipo também é típico da solução de Schwarzschild, assim como 
o decaimento do potencial efetivo, V.r(7*), igualmente em forma de lei de potência (veja 
[34]). 


Capitulo 7 


Conclusoes 


O procedimento descrito em [5] nos conduziu a uma família de soluções das equações 
de Einstein generalizadas numa brana com vácuo e constante cosmológica, A4p, não 
nula. A partir de A(r) = 1 — 2m/r — a? /r?, obtivemos apenas a métrica Schwarzschild- 
de Sitter extrema, A(r) = B(r), definida no intervalo de interesse r} < r < co, e, a 
partir de A(r) = 1 — 2m/r + a?/r?, obtivemos, numericamente, toda uma família de 
soluções do tipo Schwarzschild anti-de Sitter. Para cada C originário de (4.29), temos 
uma solução e um B(r) diferente. O nosso interesse foi direcionado para valores de C 
que tornam A(r) Æ B(r), e, para estes valores, percebemos que o comportamento assin- 
tótico, quando r cresce, desta família de soluções é o mesmo da geometria Schwarzschild 
anti-de Sitter, ou seja, quando r cresce a função B(r) tende à função A(r). Esse fato 
pôde ser ilustrado no gráfico 5.5, que mostra a diferença entre A(r) e B(r) indo a zero 
quando r cresce e no gráfico 5.6, que mostra o comportamento assintótico da coordenada 


tartaruga, 7%. 


Mostrou-se, analiticamente, que o valor da raiz da função g(r) = 4A(r) + rA,(r), rs, 
é o novo valor assumido pela singularidade física quando A(r) Æ B(r), isto é, a singu- 
laridade é transladada de r = 0 para r = rg. Numericamente foi observado este mesmo 


comportamento como pode ser visto na figura 5.4. 


O comportamento desta família de soluções, próximo ao horizonte, foi estudado 
com a obtenção da família de funções B(r) analítica. Com isso, conseguimos obter uma 
expressão analítica para a gravidade superficial, k}, em função dos parâmetros m,a 
e C. Além disso, a ordem do horizonte foi obtida com auxílio da coordenada quase 


global, u. Para C = 0, temos um horizonte de eventos de ordem 2 e para C 4 0 temos 


TT 
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um horizonte de eventos de ordem 1. 


Por fim, submetemos essas geometrias obtidas a uma perturbação escalar, ou seja, 
dado um campo escalar não massivo, ®, estudamos sua evolução nessa família de 
soluções, geometrias. O comportamento do potencial efetivo, V.y, próximo ao horizonte 
de eventos, depende, exclusivamente, do valor da constante C. Mostrou-se, analitica- 
mente, que o potencial efetivo em função da coordenada tartaruga, 7*, tende a zero 
exponencialmente quando C # 0 e tende a zero como uma lei de potência quando 
C = 0. Este resultado também foi obtido com auxílio do cálculo numérico. Toda 
família de soluções mostrou-se estável, ou seja, o campo escalar P tende a zero quando 
o tempo aumenta. Para C £ 0, o campo escalar decai exponencialmente, enquanto que 


para C = 0 o campo decai como uma lei de potência. 
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